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SUMMARY 


This  study  was  conducted  to  provide  a  database  from  which  to  draw  con¬ 
clusive  results  on  the  efficiencies  of  seismic  wave  propagation  in  natural 
terrain  and  the  resolution  and  fidelity  of  multiple  frequency  signals,  and  to 
supplement  data  for  validation  of  theoretical  models  of  seismic  wave  propaga¬ 
tion.  An  extensive  test  program  was  conducted  at  White  Sands  Missile 
Range,  New  Mexico,  using  an  electrohydraulic  vibrator,  an  impulse  loader,  and 
a  vehicle  as  sources  of  seismic  waves  over  a  5-m  to  1-km  range  and  using 
explosive  seismic  sources  over  a  1-  to  10-km  range. 

The  results  are  presented  for  discrete  frequency  vibration  tests 
(1-120  Hz),  tone  burst  tests  (1-120  Hz),  random  noise  vibration  tests,  and 
background  noise  tests  as  well  for  vehicle,  impulse,  and  explosive  tests. 
Analysis  of  data  has  been  performed  to  correlate  frequency,  amplitude,  range, 
and  other  signal  characteristics  with  model  predictions  and  to  modify  model 
coefficients  to  produce  better  predictions  for  future  tests. 

The  main  conclusions  of  this  study  are; 

a.  Seismic  signals  are  dispersed  and  selectively  attenuated  such  that 
very  low  (2-  to  3-Hz)  signals  predominate  in  the  4-  to  10-km  range, 
while  at  nearer  ranges  (less  than  4  km)  signal  frequencies  to  above 
100  Hz  are  measurable. 

b.  The  frequency  resolution  of  seismic  signals  measured  out  to  10  km  is 
below  0.5  Hz,  but  is  dependent  on  the  length  of  the  record  and  the 
stability  of  the  source. 

c.  Wind  noise  creates  problems  in  making  accurate  seismic  measurements 
over  long  ranges  (up  to  10  km)  for  wind  speeds  above  8  m/sec. 

d.  Electrical  noise  can  cause  distortions  in  measured  data  even  when 
attempts  are  made  to  remove  this  noise  with  filters. 

'  e.  The  WES  seismic  mbdel  is  shown  to  be  accurate  in  concept  for- 

defining  surface  wave  propagation  characteristics  and  for  making 
predictions  of  time  and  frequency  domain  signal  amplitudes. 

f .  Rayleigh  waves  at  low  frequency  travel  much  faster  than  at  high  fre¬ 
quency  (i.e.,  dispersion);  the  2-  to  3-Hz  component  is  controlled 
primarily  by  the  characteristics  of  deep,  high-velocity  layers  and 
the  100-  to  200-Hz  components  are  controlled  primarily  by  the  char¬ 
acteristics  of  the  slower  velocity  near  surface  layers. 

g.  Signals  measured  in  a  rock  outcrop  indicate  that  the  layers  of  soil 
above  bedrock  attenuate  the  seismic  wave  and  do  not  cause  a  focusing 
effect. 


PREFACE 


The  study  reported  herein  was  conducted  by  the  US  Army  Engineer 
Waterways  Experiment  Station  (WES)  during  Fiscal  Year  1983.  The  study  was 
funded  by  the  Ballistic  Missile  Office,  Norton  Air. Force  Base,  California, 
under  MIPR  No.  FY7653-83-00010;  the  Ballistic  Missile  Defense  System  Command, 
Huntsville,  Alabama,  under  MIPR  No.  W31RPD-33-V062;  and  the  US  Army  Corps 
of  Engineers  under  DA  Project  No.  4A762730AT42,  "Terrain  Operations  Simula¬ 
tions,"  Task  B/E3,  "Analytical  Techniques  for  the  Design  and  Application  of 
Sensors." 
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Mr.  Jerry  R.  Lundien,  Chief  of  the  Battlefield  Environment  Group  (BEG).  This 
study  was  conducted  by  Dr.  Ben  L.  Carnes  of  the  BEG.  Dr.  Carnes  and 
Mr.  Lundien  prepared  this  report. 
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to  Medium-Range  Seismic  Attenuation  Tests  Using  a  Multilayered 
Viscoelastic  Seismic  Propagation  Model,"  Technical  Report 
EL-84-9,  US  Army  Engineer  Waterways  Experiment  Station, 

Vicksburg,  Miss. 
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AN  ANALYSIS  OF  SHORT-  TO  MEDIUM-RANGE  SEISMIC  ATTENUATION  TESTS 


USING  A  MULTILAYERED  VISCOELASTIC  SEISMIC  PROPAGATION  MODEL 


PART  I:  INTRODUCTION 

Background 

1.  The  components  of  Rayleigh  waves  are  commonly  used  as  the  energy 
conveyance  whi«.h  activates  seismic  sensors.  The  amount  of  energy  generated  by 
a  target  of  military  interest  such  as  a  vehicle,  military  equipment,  or  a  spe¬ 
cific  military  operation  or  activity  is  site  dependent,  and  the  effect  of  ter¬ 
rain  on  Rayleigh  wave  propagation  is  known  to  be  substantial.  However,  little 
is  known  about  the  relative  quantitative  effects  of  specific  terrain  parame¬ 
ters  on  this  phenomenon.  In  addition,  little  is  known  about  the  propagation/ 
attenuation  of  Rayleigh  waves  over  long  ranges,  even  if  the  terrain  is  rela¬ 
tively  homogeneous  (Williams  1981).  In  an  effort  to  characterize  Rayleigh 
wave  propagation  and  attenuation  under  such  conditions,  the  US  Army  Engi¬ 
neer  Waterways  Experiment  Station  (WES)  was  requested  to  participate  in  an 
extensive  seismic  test  program.  This  program  coupled  model  calculations  with 
field  testing  to  enable  both  mathematical  and  empirical  analyses  of  wave 
propagation  and  attenuation  over  la.ge  areas. 

2.  A  mathematical  model  developed  at  WES  predicts  microseismic  signals 
in  various  terrain  materials  (Lundien  an d  Nikodem  1973).  The  model  has  been 
adapted  to  run  on  an  IBM  4331  computer  recently  installed  in  the  Environmental 
Systems  Division,  Environmental  Laboratory,  WES,  and  the  model  was  used  for 
modeling  and  predicting  the  seismic  test  site.  Predictions  are  made  rapidly 
and  efficiently  with  the  computer  model  and,  as  field  data  are  processed,  the 
model  is  tuned  to  refine  its  calculations.  When  the  model  is  refined  by  com¬ 
parison  with  measured  data,  the  model  can  be  used  to  study  the  effects  of 
other  terrain  parameters  on  the  individual  facets  of  Rayleigh  wave  generation 
and  propagation  phenomena  and  the  most  sensitive  parameters  can  be  identified. 
This  flexible  tool  can  provide  an  effective  means  of  predicting  seismic 
signals  in  are~s  where  few  data  have  been  measured. 


Objectives 


3.  The  objectives  of  this  study  were  to: 

a.  Determine  the  relative  efficiencies  with  which  seismic  signals 
propagate  from  impulsive  and  continuous  wavfe  sources  in  natural 
terrain. 

b.  Determine  the  resolution  and  fidelity  of  multiple  frequency  sig¬ 
nals  propagating  in  natural  terrain. 

c.  Supplement  an  existing  database  that  in  turn  will  be  used  to 
validate  theoretical  models  of  seismic  wave  propagation  in 
natural  terrain. 

d.  Provide  a  database  in  a  form  that  can  be  related  to  various 
military  programs  concerned  with  the  correlation  of  seismic 
information  for  detection,  classification,  and  deception 
purposes. 

Scope 

4.  In  this  study  an  extensive  amount  of  seismic  data  were  collected, 
processed,  and  compared  to  model  calculations.  Analysis  Is  limited  to  cal¬ 
culation  of  transfer  functions  and  spectral  densities  and  comparison  of  the 
same  parameters  used  in  the  model.,  Many  other  calculations  and  algorithms 
were  not  applied  because  of  limitations  in  the  scope  of  work,  but  can  be  done 
in  the  future  when  other  applications  warrant.  In  no  way  is  the  analysis  cf 
the  data  set  presented  in  this  report  complete.  '  Many  of  the  wide-ranging 
analytical  techniques  appropriate  for  use  with  these  data  are  still  under 
development . 
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PART  II:  DESCRIPTION  OF  EQUIPMENT 

Mechanical-Source  Specifications 


Vibrator 

5.  The  WES  electrohydraul ic  vibrator  is  a  portable  self-contained, 
trailer-mounted  unit  (Figure  1).  It  has  a  680-kg  inertial  mass  that  is 
excited  by  an  electronydraulic  ram.  The  system  is  capable  of  ±0.91  metric  ton 
peak  force  output  in  the  frequency  range  from  15  to  150  Hz.  A  variety  of 
forcing  functions  may  be  used,  including  steady-state  fine,  swept  sine,  tone 
burst,  and  random  excitation. 


Figure  1.  WES  electrohydrauiic  vibrator 


6.  The  vibrator  hydraulic  power  supply  is  driven  by  a  30-hp,  four- 
cylinder  gasoline  engine,  and  the  hydraulic  pump  can  provide  up  to  1.58  £/sec 
at  a  service  pressure  of  2,109  kPa.  Intercoolers  are  included  in  the 
hydraulic  system  so  that  the  vibrator  can  be  operated  continuously*  in  ambient 
temperatures  between  -4*  and  +52*  C. 

7.  In  order  to  provide  a  clean  wave  form  and  to  balance  loads  on  the 
electrohydrauiic  ram,  the  inertial  mass  is  supported  by  two  5,9-cm-diam  air 
bags.  Approximately  4.2-Pa  air  bag  pressure  is  required  to  suspend  the  inerr 
tial  mass.  The  air  bag  system  is  quite  compliant  within  the  ±0.8-ca  stroke  of 
the  ram,  so  there  is  minimal  variation  in  peak  force  output  during  shaking. 


The  electrohydraulic  ram  is  equipped  with  a  linear  variable  differential 
transformer  for  static  positioning  of  the  inertial  mass  and  for  positive  feed 
back  control  during  shaking.  The  electronic  vibrator  control  unit  compares 
the  command  signal  to  the  vibrator  with  the  feedback-  position;  signal  and 
automatically  regulates  the  ram  servovalve  to  produce  the  desired  wave  form. 
With  this  system  it  is  possible  to  produce  virtually  constant:  force  output 
from  15  to  150  Hz;  however,  the  soil  on  which  it  is  emplaced  affects  the 
actual  field  output,  as  can  be  seen  in  Figure  2.  The  vibrator  base  plate  is 
11.8  cm  in  diameter,  so  that  the  maximum  peak  dynamic  pressure  applied  to  the 
soil  at  0.91  metric  ton-force  output  is  on  the  order  of  2.0  Pa.* 


*  ,  Internal  Memorandum,  1983,  S.  S.  Cooper,  Geotechnical  Laboratory,  US  Army 
Engineer  Waterways  Experiment  Station,  Vicksburg,  Miss. 


to  increase  the  downward  velocity  of  a  falling  mass  (39-kg  solid  steel 
piston).  The  system  uses  a  small  gasoline  engine  to  power  the  vacuum  pump/air 
compressor  which  both  pressurizes  and  evacuates  air  from  the  cylinder.  In 
operating  the  impactor,  the  piston  is  first  forced  to  the  top  of  the  cylinder 
where  it  is  automatically  latched  into  position.  The  air  below  the  piston  is 
then  evacuated.  When  ready,  the  operator  pushes  a  firing  knob  which  allows 
the  piston  to  fall  1.8  m  (accelerated  by  gravity  and  air  pressure).  The 
piston  strikes  an  anvil  at  the  base  of  the  cylinder,  and  energy  is  transferred 
into  the  ground.  The  vacuum  in  the  cylinder  prevents  the  piston  from  rebound¬ 
ing  and  produces  a  single-cycle  impulsive  load.  The  impactor  is  shown  in 
Figure  3. 

T 


Figure  3.  Photograph  of  impactor 


Vehicle 

9.  The  M-35  truck  used  in  these  tests  was  a  standard  M-35  base  with  a 
shop  van  shell  mounted  on  it  (Figure  4).  The  M-35  has  tandem  axles,  each,  with 
dual  wheels  (total  of  10  wheels  on  the  ground,  counting  front)  and  weighs 
7  metric  tons  in  the  van  configuration  (Field  Manual  55-15} (Headquarters, 
Department  of  the  Army  1968).  This  vehicle  can  be  seen  in  Figure  4. 


Figure  4.  M-35  txuck  with  van  body 


Explosive-Source  Specifications 

10.  The  explosive  used  in  the  seismic  attenuation  tests  (SAT)  at  White 
Sands  Missile  Range,  New  Mexico,  was  ammonium  nitrate/fuel  oil  (ANFO) ,  which 
is  a  combination  of  prilled  ammonium  nitrate  with  about  6  percent  by  weight  of 
fuel  oil  (diesel  fuel)  as  a  sensitizer.  The  ANFO  explosion  makes  a  good 
seismic  source  with  characteristics  similar  to  those  of  TNT.  It  is  safe  and 
easy  to  use  because  of  its  relative  insensitivity  to  .temperature  and  shock. 
Explosive  charges  are  described  further  in  the  subsection  "Explosive-source 
tests,"  Part  III. 

Sensor  Specifications 

11.  The  seismic  refraction  geophones  used  in  the  test  program  were, Mark 
Products  Model  L-10,  a  standard  vertical  geophone  used  for  seismic  refraction. 
The  Model  L-10  has  a  natural  frequency  of  8.0  Hz. 

12.  The  geephones  used  in  the  seismic  tests  are  Geospace  HS-10-1  cali¬ 
brated  vertical  geophones  with  a  natural  frequency  of  1  Hz,  a  response  of 
1-200  Hz,  damped  at  70  percent  of  critical,  and  a  sensitivity  of  2.95  V/cm/ 
sec.  The  triaxial  (three  mtitualiy  perpendicular  velocities)  gedphone  was  a 
Mark  Products  Model  L-4O30  with  a  natural  frequency  of  1  Hz,  damped  70  per¬ 
cent  of  critical,  a  response  of  1-200  Hz  and  greater^  and  a  sensitivity  of 
2.35  V/cm/sec. 

13.  The  microphone  used  in  the  tests  was  a  Bruel  and  Kjaer  (8&K) 

Model  4921  outdoor  microphone  unit  with  a  frequency  range  of  10-20,000  Hz  and 
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a  dynamic  range  of  140  db.  This  unit  has  a  built-in  electrostatic  calibrator 
giving  90-db  sound  pressure  level  at  1000  Hz. 


14.  The  array  instrumentation  system  operates  from  a  commercial 

120-V  AC  power  source  that  could  not  be  provided  by  a  portable  generator  be-  • 
cause  of  the  noise  it  would  produce  in  the  array  area.  The  portable  source 
instrumentation  system  was  moved  to  each  new  source  location  and  was  posir 
tioned  with  the  same  orientation  to  the  source  at  each  location.  The  source 
instrumentation  system,  which  was  housed  in  the  van  seen  in  Figure  5a,  is 
shown  in  Figure  5b.  This  system  was  generator  powered  for  convenience.  A 
schematic  of  the  basic  analog  instrumentation  system  is  shown  in  Figure  6. 

Both  systems  described  above  were  similar  to  the  schematic,  except  the  array 
system  included  a  microphone  added  to  the  input  as  well  as  a  low-pass  filter 
between  the  amplifier  and  the  recorder.  The  lateral  geophone  system  was 
battery  powered  and  had  only  three  channels.  It  was  a  completely  self- 
contained  system  with  amplifier  and  analog  tape  recorder.  Its  geophone 
frequency  response  is  1-200  Hz;  its  natural  frequency  was  1  Hz  with  a 
sensitivity  of  2.35  V/cm/sec. 

Recorders 

15.  Signals  from  the  array  geophone  system  were  recorded  on  a  Sangamo 
Model  3500  14-trac,k  analog  tape  recorder,  recording  at  3-3/4  in.  (9.5  cm)/sec, 
FM  mode,  which  gives  a  frequency  response  of  DC  to  1250  Hz.  This  recorder 
requires  120  V  AC  power.  The  source  geophone  recorder  was  a  Racal  Model 
STORE  7DS  seven-channel  analog  tape  recorder  operated  in  FM  mode  at  3-3/4  in./ 
sec  (DC  to  1250-Hz  frequency  response).  This  recorder  was  powered  by  a  port¬ 
able  generator  operating  on,  115  V  AC. 

Amplifiers 

16.  The  array  system  employed  Ithaca  Model  456  amplifiers  with  a  band¬ 
width  of  1-100,000  Hz  and  100-db  gain  range  in  1-db  steps,  operated  from  a 
120-V  AC  power  source.  The  source  system  amplifiers  were  WES-made  units,  with 
a  bandwidth  of  DC  to  5000  Hz  and  a  continuous  gain  range  from  0  to  60  db,  and 
require  12-V  DC  power. 


Figure  5.  Photographs  of  seismic  sources  and  inscrutnentacion 
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Figure  6.  Schematic  of  typical  analog  recording  system 

Filters 

17.  The  anay  system  used  Krohn-Hite  Model  3323  filters  which  have  the 
options  for  low-pass  or  band-pass  operation.  These  filters  have  a  frequency 
range  of  0.01  to  100,000  Hz  and  an  attenuation  rate  of  24  db  per  octave  out¬ 
side  the  pass-band.  ,  The  source  geophones  were  not  filtered  since  the  signals 
were  quite  strong  at  the  source  compared  to  background  noise.  The  filters 
were,  as  a  rule  of  thumb,  set  (low-pass)  at  a  frequency  10  Hz  above  the  high¬ 
est  frequency  of  interest.  This  kept  the  filters  from  affecting  data  of  in¬ 
terest,  since  a  point  of  5  Hz  below  the  filter  setting  will  lose  approximately 

2  db  of  gain.  (Actual  filter  values  used  for  each  test  are  listed  in  Table  3.), 
Calibration 

18,  Substitution  sine-wave  calibrations  were  used  as  an  alternate  input 
to  the  amplifiers  during  system  calibration.  The  sine-wave  calibration  signal 
was  monitored  with  a  precision  voltmeter  and  frequency  counter.  The  sine-wave 
calibration  voltage  was  compared  to  the  known  geophone  sensitivities  to  cal¬ 
culate  the  recording  system  sensitivity  in  velocity  units.  A  time  standard. 


IRIG  B,  was  recorded  on  both  tape  recorders  for  time  of  day  and  common  timing 
on  both  tape  recorders. 

19.  Calibration  of  all  geophone  channels  was  accomplished  with  a  sub¬ 
stitution  sine-wave  signal  input  to  the  amplifier.  The  sine  wave  was  measured 
with  an  accurate  voltmeter  and  proportioned  to  a  velocity  value  by  knowing  the 
sensitivity  of  each  geophone.  The  amplifier  output  was  connected  to  the 
Krohn-Hite  filters,  which  were  connected  to  the  tape  recorder.  A  30-Hz  sine 
wave  was  selected  as  the  calibration  frequency. 

20.  During  calibration,,  the  smallest  setting  was  selected  on  each 
amplifier,  approximately  15  db,  with  the  remote  gain  control.  An  input  signal 
of  125  mV  root  mean  square  (RMS)  amplitude  was  applied  to  each  amplifier 
simultaneously,  through  the  Krohn-Hitc  filter,  set  at  140-Hz  low  pass  and  on 
to  the  tape  recorder.  The  125-mV  RMS  signal  amplified  with  a  net  gain  of 

+12  db  gives  0.5-V  RMS  on  the  tape  recorder.  One  minute  of  30-Hz  signal  was 
recorded.  On  several  calibrations,  a  sine-wave  sweep  from  1  Hz  to  200  Hz  was 
also  recorded. 

21.  Each  geophone  channel  was  filtered  before  being  recorded  by  the 
Krohn-Hite  filters.  Each  unit  held  two  separate  filter  circuits,  and  two  cas¬ 
caded  low-pass  filters  were  used  on  each  geophone  channel.  This  gave  an 
attenuation  above  the  corner  frequency  of  48  db  per  octave.  The  corner  fre¬ 
quency  was  adjusted  for  each  test,  with  a  rule  of  thumb  that  the  corner  fre- 

c-  icy  be  above  the  vibrator  driver  frequency  by  10  Hz. 

22.  The  amplifier  gain  was  adjusted  for  each  test  such  that  ''he  input 

i  » 

to  the  tape  recorder  during  a  test  would  be  less  than  1  V  RMS.  All  geophone 
amplifiers  were  set  to  provide  the  1-V  RMS  maximum  input  for  each  .test  with 
the  use  of  the  remote  gain  control.  .. 

23.  The  B&K  microphone  channel  was  calibrated  using  a  built-in  signal 
and  electrostatic  calibrator  mounted  directly  on  the  microphone.  The  calibra-, 
tion  frequency  was  1000  Hz,  and  the  signal  was  90-db  sound  pressure  level.  No 
filtering  was  used  on  the  B&K  microphone  signal. 

Spectrum  analyzer 

24.  The  spectrum  analyzer  was  used  with  the  array  instrument 3 ♦'ion  sys¬ 
tem  to  make  field  determinations  as  to  what  and  how  background  noise  data 
could  be  filtered  from  actual  test  signals.  It  was  also  used  to  determine 
during  tasting  if  various  input  frequencies  were  present  in  the  signals,  and 
what  frequencies  were  uominant  in  impulse  and  explosive  tests.  The  spectrum 


analyzer  used  was  a  Hewlett-Packard  3582,  which  has  a  frequency  response  of  DQ 
to  20,000  Hz  and  a  dynamic  range  of  70  db. 


PART  III:  DESCRIPTION  OF  TEST  PROGRAM 


25.  Upon  identification  of  objectives  and  scope  by  the  sponsor,  an 
extensive  test  plan  was  devised  which  included  proper  and  sufficient  tests  to 
accomplish  the  objectives;  the  test  plan  is  presented  in  its:  entirety  in  Ap¬ 
pendix  A.  The  first  activity  was  to  select  a  site  and  conduct  a  preliminary 
site  characterization.  The  basic  site  requirements  were:  reasonable  uniform¬ 
ity  in  subsurface  conditions  over  10  km,  flat  to  rolling  topography,  depth  to 
water  table  at  least  42  m,  and  seismic  profile  beginning  at  427  m/sec  down  to 
6  m,  900-1370  m/sec  between  6  m  and  42  m,  and  ±1500  m/sec  below  42  m.  The 
site  requirements  are  described  further  in  Appendixes  A  and  B.  After  consid¬ 
eration  of  numerous  sites  within  the  continental  United  States,  two  potential 
sites  were  selected,  based  on  available  data.  One  was  near  the  Oscuro  gate  in 
the  northeast  area  of  the  US  Army  White  Sands  Missile  Range  (WSMR) ,  New 
Mexico.  The  second  was  near  the  Trinity  site  in  the  northwest  area  of  WSMR. 
The  Oscuro  site  was  the  most  desirable  because  of  low  background  noise.  Dur¬ 
ing  preliminary  site  characterization  tests,  however,  a  layer  of  rock  with 
compression  wave  velocity  on  the  order  of  3000  m/sec  was  found  within  10  m  of 
the  surface,  and  the  site  was  deemed  unusable  for  the  tests  because  the  shal¬ 
low  rock  layer  did  npt  meet  the  requirements.  The  second  site  was  verified  as 
being  usable,  with  several  seismic  refraction  lines  (included  in  Appendix  B) 
in  correlation  with  previous  data  taken  in  the  same  area.*  The  site,  although 
seismically  compatible  with  test  requirements  (see  Appendix  3).,  did  pose  some 
operational  problems  with  periodic  increases  in  background  noise  from  con¬ 
struction  activity  In  the  area  for  DIRECT  COURSE,  a  major  weapons  test. 

26.  The  second  preliminary  activity  for  this  test  program  was  to  cali¬ 
brate  and  check  the  seismic  sources  at  the  WES.  The  vibrator  and  the  impactor 
were  both  tested  at  the  WF.S  in  preliminary  tests  to  determine  operating  char¬ 
acteristics  in  the  planned  modes  of  use.  The  vibrator  was  tested  over  the 
expected  frequency  range;  these  preliminary  data  are  shown  in  Figure  2.  The 


*  Internal  Memorandum,  "Operation  DICETHROW,",  1975,  S.  S.  Cooper,  Geotechni¬ 
cal  Laboratory,  US  Army  Engineer  Waterways  Experiment  Station,  Vicksburg, 
Miss. 


impactor  was  also  tested  and  showed  a  clean  impulse  load.  Preliminary  results 
were  used  in  developing  plans  and  procedures  for  extracting  needed  data  from 
the  field  tests. 

Test  Site  ' 

27.  The  SAT  program  was  conducted  near  the  Trinity  site  at  the  WSMR. 

The  tf.st  course  layout  began  near  WSMR  Route  13  and  went  north,  passing  just 
west  of  Trinity  and  ending  near  Beck  (see  site  location  map.  Figure  7).  The' 
test  line  passed  about  1.0  km  west  of  the  Trinity  center.  Calculations  from 
weapons  test  data  reveal  that  the  Trinity  event,  a  19-kiloton  nuclear  weapon 
detonated  on  a  30- m  tower  (Rooke,  Carnes,  and  Davis  1974),  would  produce  per¬ 
manent  displacement  or  damage  to  the  substrate  out  to  a  maximum  distance  of 
300  m  (Carnes  1974).  Thus,  the  test  line  was  located  a  sufficient  distance  to 
be  in  undisturbed  material. 

28.  The  general  surface  conditions  at  the  SAT  site  were  relatively  flat 
to  rolling,  covered  by  sparse  vegetation  consisting  of  creosote  brush,  desert 
grass,  cactus,  and  yucca  which  varied  from  10  to  35  percent  (average)  over  the 
test  site.  The  site  soil  was  a  silty  sand  near  the  surface  that  gradually 
graded  coarser  with  depth.  It  was  difficult  to  travel  cross-country  over  the 
site  because  the  sand  was  loose  in  some  areas  and  accumulated  around  the  veg¬ 
etation  in  others.  Photographs  of  the  site  are  shown  in  Figure  8.  The  test 
site  lies  in  the  Jornada  del  Muerto  Basin  of  the  Mexican  Highland  Section  of 
the  basin  and  Range  Province,  which  is  characterized  by  a  series  of  fault 
blocks  forming  asymmetric  ridges  or  mountains  and  broad  intervening  basins, 
this  basin  has  an  overburden  of  alluvial  material  containing  sand  and  rock 
fragments,  generally  grading  coarser  with  depth.* 

Site  Documentation 

Seismic  refraction  tests 

29.  The  seismic  refract'on  data  from  both  the  preliminary  and  the  site 
characterization  tests  were  obtained  using  a  portable  24-channel  seismograph. 

*  Internal  Memorandum,  1983,  S.  S.  Cooper,  Geotechnical  Laboratory, 

US  Army  Engineer  Waterways  Experiment  Station,  Vicksburg,  Miss. 
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a.  View  of  array,  looking  east,  showing  B&K  microphone 


b.  View  of  10-km  location  looking  north 


Figure  8.  Photographs  of  site 
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This  seismograph  produces  a  permanent  record  on  oscillograph  paper.  A  total 
of  23  vertical,  velocity-type  geophones  were  used  in  the  conduct  of  the  tests. 
The  geophones  were  emplaced  in  a  straight  line  along  the  surface  of  the  ground 
at  50-m  intervals.  The  seismic  source  was  a  2-kg  explosive  charge  detonated 
50  m  from  one  end  geophone  and  repeated  at  the  opposite  end.  This  resulted  in 
two  seismic  traverses  (forward  and'  reverse),  366  m  in  length,  and  delineation 
of  substrate  do^n  to  a  minimum  depth  of  90  m,  in  accordance  with  Engineer 
Manual  (EM)  111 0— 1  — 1802  (Headquarters,  Department  of  the  Army  1979). 

30.  Data  obtained  from  the  surface  refraction  seismic  tests  consist  of 
the  time  required  for  a  compression  wave  to  travel  from  a  seismic  source  to 
points  of  measurement.  The  data  are  plotted  as  time  versus  distance  from 
which  compression  wave  velocities  and  depths  to  interfaces  can  be  computed  by 
using  conventional  seismology  equations  (EM  1110-1-1802).  These  data  and 
techniques  are  discussed  in  detail  in  Appendix  B.  Shear  wave  velocities  were 
not  determined  for  these  tests.  Since  the  refraction  test  lines  were  too 
short  to  determine  shear  wave  velocities  from  explosions  while  using  vertical 
geophones,  empirical  relationships  derived  from  previous  test  programs  were 
used  to  estimate  the  shear  wave  velocities  for  input  to  the  theoretical  model. 
Environmental  data 

31.  Quantitative  environmental  data  describing  the  surface  and  sub¬ 
strate  conditions  included  soil  wet  density,  soil  moisture  content,  ground 
cone  index,  soil  type,  soil  grain  size,  distribution,  surface  topography, 
vegetation  type  and  density,  surface  load-response  characteristics,  and 
meteorological  conditions.  The  data  listed  below  for  the  site  are  summarized 
in  Table  l. 

32.  Soils,  Soil  samples  representative  of  the  surface  were. obtained 
for  laboratory  analysis  by  digging  to  a  30-cm  depth.  These  samples  were 
tested  in  the  WES  soil  laboratory  for  moisture  content  and  grain  size 
distribution. 

33.  Cone  index.  Cone  index  is  a  measure  of  soil  strength  measured  with 
a  1.25-sq-cm  cone  and  expressed  herein  in  kilooascals.  It  is  determined  by 
measuring  the  force  needed  to  move  a  cone  penetrometer  vertically  into  the 
ground  at  a  constant  rate  (the  maximum  reading  available  on  the  hand  cone 
penetrometer  is  750).  Cone  index  data  were  taken  at  the  array  and  at  the 
0.1-km,  4.0-km,  and  10.0-km  test  locations,  near  the  sources  and  the  vehicle 


34.  Topography.  The  topography  of  the  site  is  shown  on  US  Geologi¬ 
cal  Survey  topographic  map  no.  N13-/ ,  (Tularosa) ,  Series  V502,  scale  1:250,000. 
However,  for  purposes  of  this  program,  the  most  important  topographic  data 
come  from  an  engineering  survey  along  the  test  line.  This  was  accomplished 
using  a  laser-type  electronic  distance  measuring  (EDM)  system  that  measures 
both  horizontal  and  vertical  distance  between  points.  In  addition  to  the  site 
profile,  profiles  of  each  vertical  track  were  measured  upon  completion  of  the 
tests.  It  should  be  noted  that  these  track  profiles  represent  an  instantan¬ 
eous  example  of  the  degree  of  roughness  that  prevails  on  a  particular  vehicle 
run. 

35.  Meteorological  cond' cions.  A  meteorological  station  was  estab^- 
lished  at  the  site  for  the  duration  of  the  tests.  Measurements  of  air  and 
ground  temperature,  wind  speed  and  direction,  solar  insolation,  relative 
humidity,  and  precipitation  were  made  continuously  and  recorded  automatically 
every  30  min.  The  data  were  recorded  on  a  cassette  and  then  transferred  to  a 
computer  for  reduction  and  processing.. 

Seismic  Wave  Attenuation  Tests 

Mechanical-source  tests 

36.  Seismic  wave  attenuation  tests  were  conducted  using  several  mechan¬ 
ical  seismic  sources,  including  a  vibrator,  an  impact  source,  and  an  M-35 
truck.  The  tests  were  conducted  in  the  layout  shown  in  Figure  9.  An  array  of 
calibrated  sensors,  including  one  three-dimensional  geophone,  five  vertical 
geophones,  and  one  outdoor  microphone  uric,  were  emplaced  ft  the  beginning  , 
point  of  the  test  line  to  comprise  the  array  sensor  system.  The  three  sources 
were  employed  individually  at  0.1  km,  0.5  km,  and  1.0  km,  but  no  data  we'c 
collected  at  ranges  longer  than  1  km  with  these  sources  since  measurable 
signals  were  at  the  level  of  background  noise  or  less  for  ranges  of  2  km  and 
beyond.  Emplaced  at  each  source  location  were  four  vert'1  sal  geophones  and  one 
outdoor  microphone  unit  to  measure  close-in  attenuation,  signal  variation 
between  locations,  and  background  noise.  In  addition,  a  load  cell  was  used  to 
measure  force  applied  to  the  ground  by  the  vibrator,  and  an  accelerometer  was 
used  to  measure,  impactor  acceleration,  from  which  ground  loading  was 
calculated. 
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37.  Vibration  tests.  The  vibrator  was  run  at  discrete  frequencies 
between  1  and  120  Hz.  as  listed  in  Table  2,  for  periods  of  approximately 

60  sec  while  data  were  recorded  at  the  source  and  at  the  array.  The  vibrator 
was  also  run  at  a  10-percent  offset  frequency  (i.e.,  for  nominal  frequency, 
f=50  Hz,  offset  frequency  f=55  Hz)  for  five  separate  frequencies.  In 
addition,  tone  bursts  were  run  at  each  discrete  frequency  with  the  vibrator  on 
for  8  to  32  cycles  and  off  for  8  to  128  cycles  (to  give  a  reasonable  cycle  of 
signal  and  quiet)  and  repeated  for  about  60  sec.  The  last  tests  were  run  with 
the  vibrator  set  with  random  noise  input  for  approximately  60  sec  and  then 
repeated.  All  tests  conducted  with  the  vibrator  are  listed  in  Table  2. 

38.  Impulse  tests.  The  impactor  was  placed  into  position  and  released 
five  times  to  seat  the  plate  and  establish  proper  gains  for  the  instrumenta¬ 
tion  systems.  It  was  then  released  five  times  for  recording  of  the  seismic 
signal  with  a  definite  break  (60  sec)  between  impacts.  The  impulse  tests 
conducted  are  also  listed  in  Table  2. 

39.  Vehicle  tests.  The  vehicle  tests  were  conducted  with  an  Army  2.5- 
ton  truck,  M-35,  as  the  source  of  seismic  and  acoustic  signals.  The  M-35 
truck  provides  a  complex  forcing  function  composed  of  both  impulsive  and 


multiple  frequency  components  (Engdahl  and  West  1974).  The  terrain  was  so 
rough  that  the  vehicle  could  safely  be  driven  to  speeds  of  only  16  km/hr 
(greater  speeds  produced  excessive  driver  motion) .  The  vehicle  was  driven  at 


both  8  km/hr  and  16  km/hr  perpendicular  to  the  test  line  for  a  distance  of 
about  50  m  on  each  side  of  the  ling,  allowing  approximately'  5m  for  accelera¬ 
tion  aud  deceleration  at  either  end.  The  vehicle  was  used  unloaded,  and  since 
the  path  was  rough,  the  vehicle  motion  and  vibration  at  both  speeds  produced  a 
forcing  function  for  seismic  wave  excitation.  This  vehicle  response  can  be 
used  to  relate  response  of  both  smaller  and  larger  vehicles  by  relating  the 
forcing  function  and  measured  signal  with  model  calculations  for  the  larger  or 
smaller  forcing  functions. 

Explosive-source  tests 

40.  Since  the  nondestructive  attenuation  tests  did  not  produce  data 
beyond' 1.0  km,  explosive  sources  were  used  to  extend  th *  data  to  longer 
ranges.  Explosive-source  tests  were  run  at  ranges  of  1,  2 ,  4,  8,  and  10  km 
(Table  3) .  The  test  at  1  km  was  run  to  obtain  data  to  relate  the  explosive- 
source  response  with  that  of, the  mechanical  sources.  Preliminary  tests  were 
designed  to  define  a  specific  explosive  charge  that  could  propagate  a 


measurable  signal  to  the  array  from  10  km.  The  first  explosive  shot,  which 
weighed  12.25  kg,  produced  a  good  signal  at  the  array  (this  shot  size  was 
selected  on  the  basis  of  pretest  predictions  during  the  planning  cycle). 

Thus,  the  same  size  charge  was  used  at  each  site  thereafter,  except  the  last 
shot  (number  136),  where  excess  explosives  weighing  63.5  kg  were  detonated 
when  the  program  was  completed.  (This  also  provided  a  check  of  seismic  signal 
for  two  different  strength  energy  sources  at  one  range.)  The  explosive  used 
in  all  the  tests  was  an  ANFO  mixture  buried  3  m  for  good  energy  coupling. 

41.  For  the  explosive  tests,  the  source  geophones  were  emplaced  in  a 
different  pattern,  since  signal  levels  were  expected  to  be  much  higher  than 
with  the  mechanical-source  tests.  Two  vertical  geophones  were  emplaced  300  m 
from  the  shot,  and  one  vertical  geophone  was  emplaced  100  m  from  the  shot 
(except  for  shot  number  136,  for  which  it  was  emplaced  at  200  ni) .  In  all  but 
the  8-km  shot  (number  135),  the  geophones  were  emplaced  between  the  shot  loca¬ 
tion  and  the  array.  However,  in  this  case  the  geophones  were  left  in  the  same 
emplacement  as  for  the  10-km  shoes  (numbers  134  and  136) ,' leaving  the  geo¬ 
phones  located  1700  m  and  1900  m  from  the  8-km  shot. 

42.  Since  the  explosive  tests  were  expected  to  propagate  a  measurable 
signal  to  the  array  over  long  distances,  an  add-on  test  was  designed  to  mea¬ 
sure  the  signal  in  an  outcrop  of  the  substrate.  One  triaxial  geophone  was 
emplaced  in  an  outcrop  lying  about  1 1  km  to  the  east  of  the  10-km  station. 

This  lateral  system  was  used  to  determine  signal  characteristics  at  a  high- 
contrast  interface.  The  geophone  was  monitored  with  a  battery-powered  seismo¬ 
graph  that  produced  oscillograph  records  of  the  measured  signals. 


PART  IV:  RESULTS 


Physical  Site  Data  Obtained 


Subsurface  properties 

43.  The  rest  site  overburden  typically  consists  of  fine  material  inter- 
bedded  with  occasional  layers  of  coarse  sand.  The  material  grades  coarser 
with  depth  and  can  be  expected  to  have  some  patches  of  caliche  and  conglomer¬ 
ate.  However,  the  seismic  refraction  tests  showed  a  reasonable  uniformity  in 
subsurface  conditions  over  the  iO-km  test  line.  The  average  seismic  P-wave 
velocities  are  listed  in  Table  1.  An  approximate  cross  section  of  subsurface 
properties  is  shown  in  Figure  10.  It  should  be  noted  that  this  cross  section 
is  a  generalized  view  of  the  site  based  on  the  seismic  refraction  data  and 
available  information,  coupled  with  the  experience  of  the  seismologist. 


Figure  10.  Seismic  profile  of  SAT  site 
Surface  soil  properties 

44.  The  surface  soil  is  a  silty  sand,  classified  SM  by  the  Unified  Soil 
Classification  System  (Lambe  and  Whitman  196$).  It  had  a  moisture  content  of 
about  (>  percent  on  the  surface  and  about  8.8  percent  at  a  depth  of  30  cm.  The 


soil  had  a  cone  index  of  100  to  220  on  the  surface  and  increased  to  approxi¬ 
mately  300  at  a  depth  of  37.5  cm.  All  measured  properties  are  listed  in 
Table  1 . 

Topography 

45.  The  site  topography  is  summarized  in  Figure  11,  which  shows  a  maxi¬ 
mum  site  elevation  difference  of  about  32  m,  gradually  increasing  in  elevation 


from  the  array  to  the  10-km  station.  There  were  a  few  localized  anomalies, 
usually  in  a  wash  location. 

Meteorological  conditions 

46.  Air  temperature  and  wind  speed  taken  at  the  meteorological  station 
for  each  test  are  listed  in  the  last  two  columns  of  Tables  4-6.  Wind  veloci¬ 
ties  on  several  days  were  qpite  strong.  These  conditions  are  contrasted  with 
those  on  calmer  days  in  Figure  12.  Similarly,  ground  and  air  temperatures 
varied  considerably  over  the  test  period,  as  shown  in  Figure  13.  Although 
some  light  precipitation  occurred  during  the  preparation  period,  no  sig¬ 
nificant  precipitation  was  experienced  during  the  tests.  Ground  temperatures 
remained  above  freezing  during  the  test  period.  Measured  weather  data  are 
summarized  in  Table  7. 


TIME,  HOURS 

Figure  12.  Diurnal  wind  speed  for  selected  days 

*  i 

I 

Seismic  Wave  Attenuation  Data 

*  Mechanical  sources 

j  47.  Vibration  data.  The  result  of  the  vibration  tests  was  a  great  vol¬ 

ume  of  data  involving  about  1500  channels  of  analog  records.  These  data  are 

;  presented  in  Plates  1-84,  as  both  time  and  frequency  domain  plots;  peak  values 

!  are  listed  in  Table  4.  The  discrete  frequency  vibration  tests  are  presented 

in  Plates  1-43.  Each  plate  presents  the  results  of  one  test  with  the  input 
load  shown  at  the  top;  next  the  signal  measured  at  50  m;  and,  continuing  down 
the  page,  the  vertical,  radial,  and  transverse  signals  measured  by  the 
triaxial  geophone  at  the  main  array.  Plots  of  time  and  frequency  domain  for 
each  signal-  are  presented  side-by-side  on  the  plates  for  a  quick  look  at  the 

*  frequency  content  of  the  signal. 

48.  The  signals  show  that  measured  frequency  is  within  the  margin  of 
measurement  error  of  the  nominal  input  frequency.  In  the  10-percent  offset 
frequency  tests,  measured  frequency  is  distinguishable  from  the  base  fre- 
quency,  even  down  to  3.3  rtz  and  3.0  Hz.  However,  because  of  the  digitizing 
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sampling  rate  cf  1000  points  per  second  and  a  total  2048  points  sampled,  there 
is  a  minimum  frequency  resolution  of  0.488  Hz  in  the  "fast"  Fourier  Trans¬ 
form  (FFT).  In  this  case  the  FFT  reads  2.928  Hz  for  the  3.0-Hz  signal  and 
3.416  Hz  for  the  3.3-Hz  signal,  which  are  reasonable  values.  However,  to 
check  frequency  differences  below  0.25  Hz,  a  different  sampling  rate  would  be 
required.  There  appears  to  be  little  difference  in  amplitude  between  ,the 
discrete  and  offset  frequency. 

49.  The  tone  burst  test  data  presented  in  Plates  44-77  are  in  a  form 
similar  to  tlie  discrete  frequency  tests  with  time  and  frequency  domain  plots 
of  input  lead,  50-m  source  geophone,  and  triaxial  array  geophone  records.  The 
tone  bursts  give  an  opportunity  to  describe  the  accuracy  of  the  measurements 
at  each  recording  station  by  comparing  amplitude  and  frequency  results  with 
corresponding  discrete  frequency  test  results.  The  measured  amplitudes  on  the 
tone  burst  tests  compare  favorably  with  the  discrete  frequency  results.  There 
are  some  differences  in  the  peak  velocities  as  measured  in  the  FFT  plots, 
probably  caused  by  a  "rounding"  of  the  input  load  and  by  the  introduction  of 
harmonics  by  the  vibrator  in  a  tone  burst  mode.  The  difference  does  not  nec¬ 
essarily  seem  to  be  dependent  on  frequency  or  the  number  of  cycles,  but  rather 
on  whether  the  sinusoidal  input  is  "clean."  There  are  some  other  inherent 
p-oblems  in  running  the  vibrator  in  a  tone  burst  mode.  At  frequencies  below 

5  Hz,  the  vibrator  turns  on  and  off  with  a  much  larger  load  than  during  the 
"on"  cycles  of  the  tone  burst.  The  result  is  a  "framed"  tone  burst,  as  seen 
in  Plates  55,  56,  and  57.  It  should  be. noted  that  the  FFT  was  selectively 
run,  omitting  these  peak  load  sections  of  the  signal. 

50.  The  frequency  sweep  tests  are  presented  in  Plates  78-81.  The  time 
domain  plot  presented  in  Plate  78  for  the  1-km  sweep  test  shows  the  response 
of  the  site  at  the  5-,  20-,  50-,  100-,  and  1000-m  ranges  to  the  vibrator  fre¬ 
quency  sweep.  The  time  domain  plot  for  the  1-km  site  is  presented  as  repre¬ 
sentative  of  the  sweep  tests  since  the  tests,  which  took  about  4  min  to  com¬ 
plete,  are  difficult  to  display.  The  FFT  plots  of  Plates  79-81  (for  all 
locations)  show  amplitude  versus  frequency.  The  plot  in  Figure  14  is  a  com¬ 
posite  of  the  average  peak  FFT  for  each  range  (taken  from  Plates  80  and  81), 
showing  the  variation  of  signal  with  both  frequency  and  range. 

51.  Results  of  the  random  noise  tests  are  presented  in  Plates  82-84. 

The  data  from  the  source  geophones  are  good.  However,  the  data  from  the  afray 
sensors  are  affected  by  the  60-Hz  electrical  noise,  and  even  with  a  40-Hz 
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low-pass  filter,  there  is  a  large  60-Hz  component  in  the  signal.  The  data 
from  the  500-m  and  1-km  ranges  probably  are  not  affected  by  the  40-Hz  filter, 
but  the  iGO-m  data  are  affected  somewhat  since  some  40-.  to  120-Hz  data  are 
present  in  the  signal  at  100  tn.  Thus,  the  100-m  source  geophone  data  are 
substituted  for  the  100-m  vertical  data  from  the  array. 

52.  Impulse  data.  The  impulse  test  data  are  presented  in  Plates  85-87 
and  peak  values  are  listed  in  Table  5.  Each  plate  represents  one  test,  with 
the  impulsive  acceleration  converted  to  an  approximation  of  load,  at  the  top;  , 
the  signal  measured  at  50  m,  next;  and  then  the  vertical,  radial,  and  trans¬ 
verse  signals  measured  by  the  triaxial  geophone  at  the  main  array.  Both  time 
and  frequency  domain  plots  are  shown  for  each  signal.  These  data  show  the 
effect  of  the  site  on  frequency  content  of  the  signal,  as  the  dominant  fre¬ 
quency  decreases  with  increasing  range. 

53.  Vehicle  data.  The  vehicle  data  are  presented  in  Plates  88-93. 

Each  plate  represents  one  test,  with  the  5-m  geophone  data  at  the  top;  fol-' 
lowed  by  the  50-m  geophone  data;  and  then  the  vertical,  radial,  and  transverse 
components  of  the  triaxial  geophone  at  the  main  array.  Both  slow  and  fast 
vehicle  speed  tests  are  presented,  and  both  time  and  frequency  domain  plots 
are  shown  for  each  signal. 

Explosive  sources 

54.  Data  from  the  explosive-source  seismic  attenuation  tests  are  pre¬ 
sented  in  Plates  94-99  and  peak  values  are  listed  in  Table  6.  Each  plate 
represents  one  test  with  the  close-in  data  for  100  m  at  the  top;  followed  by 
the  300-m  data;  and  then  the  vertical,  radial,  and  transverse  components  of 
data  measured  by  the  triaxial  geophone  at  the  main  array.  Time  and  frequency 
domain  plots  of  each  signal  are  shown  side-by-side.  The  data'  are  generally 
good,  although  there  is  some  wind  noise  interference  at  the  main  array  on  some 
of  the  tests,  as  can  be  seen  in  Plates  97  and  98  (Tests  135  and  134).  The 
8-km  test  data,  in  particular,  illustrate  how  wind  noise  obscures  the  signal 
(signal  is  the  same  before  and  during  the  active  test  period),  and  since  fre-. 
quency  contents  are  similar  between  the  two  types  of  signal,  the  data  cannot 
be  easily  filtered  to  remove  the  noise. 

55.  The  data  for  the  lateral  geophone  system  are  limited  to  the  explo¬ 
sive  tests.  The  system  was  emplaced  in  a  position  about  11  km  perpendicular 
to  the  test  line  (see  site  location  map.  Figure  7)  and  was  calibrated  at  the 
same  levels  as  the  array  geophones  for  the  10-km  tests.  Timing  was  difficult 
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because  of  the  limited  communication  afforded  by  portable  radios  with  short 
ranges . 

56.  The  data  were  limited  to  those  recorded  on  a  strip  chart  recorder, 
with  less  fidelity  than  data  recorded  on  a  tape.  However,  they  address  the 
primary  question,  whether  any  unusual  seismic  signal  was  loaded  into  bedrock 
by  the  explosive.  The  triaxial  geophon®  was  see  into  a  rock  outcrop  in  this 
lateral  position  co  give,  an  indication  of  any  unusual  seismic  signals.  Since 
the  geophone  was  calibrated  for  the  same  signal  at  the  array,  a  large  signal 
would  indicate  some  unusual  loading.  However,  no  signal  received  at  the 
lateral  system  was  as  large  as  the  array,  and  most  signals  were  significantly 
smaller.  As  can  be  seen  in  Figure  15,  the  signal  produced  by  Test  132 


Figure  15.  Envelope  records  from  lateral  geophone  system  showing 
maximum  and  minimum  values  for  each  trace  vTest  131) 
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produced  a  maximum  amplitude  of  about  75  x  10  ^  cm/sec  while  the  signal  for 
the  explosive  tests  received  at  the  array  did  not  fall  below  700  x  10  ^  cm/sec. 
Thus,  the  bedrock  was  deep  enough  at  the  site  to  be  of  little  significance. 


PART  V:  ANALYSIS  AND  DISCUSSION 


Seismic  Waves 

57.  The  generation,  propagation,  and  sensing  of  seismic  energy  from 
targets  of  military  interest  consist  of  a  sequence  of  interrelated  physical 
phenomena  (see  Figure  16) .  Anything  that  moves  on  the  surface  of  the  ground 
applies  a  force  to  the  substrate  which  can  be  sensed  at  some  distance  from 
that  source.  If  the  force  is  large  enough  to  stress  the  substrate  to  a 
dynamic  deformation,  the  energy  is  carried  from  that  point  by  seismic  waves. 
The  amplitude  and  type  (to  some  degree)  of  seismic  wave  generated  depend  on 
the  force  history  applied  by  the  target.  Variations  in  the  target  which  pro¬ 
duce  changes  in  the  force  history  are  mainly  mass  and  rate  of  force  applica¬ 
tion.  A  large  mass  will  cause  greater  forces  than  a  small  mass  and  will 
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Figure  16.  Properties  of  seismic  attenuation 


result  in  a  different  force  history.  In  addition,  a  given  force  applied 
quickly  will  produce  a  different  force  history  than  force  applied  slowly. 
Variations  in  the  nature  of  the  substrate  may  also  affect  the  force  history, 
since  that  history  is  a  record  of  the  relation  between  force  and  time  as 
measured  at  a  specific  point. 

58.  There  are  three  major  seismic  wave  modes,  all  of  which  move  radi¬ 
ally  outward  from  a  point  of  substrate  deformation  (i.e.  the  point  at  which 
stress  is  applied):  (a)  compression  waves  (P-waves),  in  which  the  principal 
particle  motion  is  along  a  radial;  (b)  shear  waves  (S-waves),  in  which  the 
principal  particle  motion  is  at  right  angles  to  a  radial;  and  (c)  surface 
waves  (Rayleigh  waves),  in  which  the  principal  particle  motion  is  elliptically 
retrograde  in  a  plane  chiefly  perpendicular  to  the  surface  of  the  propagating 
medium.  Compression  waves  move  in  all  directions,  so  that  the  advancing  wave 
front  is  approximately  a  hemisphere.  However,  they  do  not  necessarily  move  at 
the  same  velocity  as  the  S-wave,  even  in  the  same  material.  Unlike  P-waves 
and  S-waves,  Rayleigh  waves  generally  move  at  shallower  depths  and  parallel  to 
the  surface,  so  that  at  a  point  a  few  metres  from  the  target,  the  advancing 
wave  front  is  cylindrical  with  the  vertical  axis  of  the  cylinder  at  the  tar¬ 
get.  For  all  practical  purposes,  Rayleigh  waves  (assuming  the  very  low  ampli¬ 
tudes  that  are  of  interest  in  this  study)  affect  a  substrate  depth  equal  to 
about  one-third  wavelength.  Amplitudes  of  Rayleigh  waves  are  usually  at  a 
maximum  at  the  surface  and  decrease  progressively  with  depth.  The  propagation 
velocities  are  not  the  same  as  those  of  either  P-waves  or  S-waves. 

59.  Generally,  boith  P-waves  and  S-waves  are  diffracted  in  the  direction 
of  greater  substrate  densities,  and  since  nearly  all  natural  substrates  in¬ 
crease  in  density  with  depth,  these  waves  tend  to  diffract  downward.  The 
result  is  that  they  appear  to  attenuate  very  rapidly  at  the  surface.  They 
may,  however,  reflect  from  a  subsurface  discontinuity  and  reappear  at  the  sur¬ 
face  at  some  distance  ftom  the  target.  Thus,  there  could  be  zones  in  which 
P-waves  and  S-waves  are  difficult  to  detect  at  the  surface.  Conversely,  the 
Rayleigh  wave  generally  propagates  along  the  surface;  thus,  theoretically, 
there  can  be  no  discontinuities  in  the  area  over  which  the  signal  can  be 
detected.  This  knowledge,  in  addition  to  the  fact  that  about  two-thirds  of 
the  energy,  at  the  source  is  carried  away  in  the  Rayleigh  wave,  is  the  princi¬ 
pal  reason  for  choosing  the  Rayleigh  wave  as  the  mode  upon  which  to  base  this 
study. . 
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60.  The  "ideal"  or  theoretical  situation  for  the  propagation  of 
Rayleigh  waves  would  consist  of  a  completely  homogeneous  elastic  half-space. 
Anything  that  departs  from  these  ideal  conditions  serves  to  attenuate  or  dis¬ 
perse  the  wave  form  more  rapidly  than  .would  be  the  case  by  purely  geometric 
attenuation  (i.e.,  the  attenuation  resulting'  from  the  same  amount  of  energy 
being  applied  over  a  longer  wave  front  as  the  yave  moves  outward  away  from  the 
source) .  There  are  several  basic  types  of  nonhomoger  ies  commonly  found  in 
nature.  The  dramatic  effects  of  nonhomogeneities  on  *ave  propagation  stem 
primarily  from  the  fact  that  each  type  of  material  tends  to  act  as  a  specific  . 
medium  for  a  specific  suite  of  wave  frequencies.  That  is,  a  specific  type  of 
material  tends  to  propagate  certain  frequencies  more  efficiently  than  othpre. 
The  result  is  that  substrate  materials  act  as  selective  filters.  For  example, 
most  targets  generate  a  seismic  wave  train  (or  signal)  consisting  of  a  complex 
of,  frequencies  (or  wavelengths)  ranging  from  very  low  (i.e,,  very  long  waves) 
to  very  high  (i.e.,  very  short  waves).  Generally,  the  signal  contains  a  broad 
spectrum  of  frequencies  as  it  emerges  from  the  source.  However,  as  the  wave 
train  moves  radially  away  from  the  source,  two  things  happen  to  it: 

a.  Since  each  frequency  tends  to  propagate  at  a  slightly  different 
speed,  the  wave  train  tends  to  separate  into  sections  (dis¬ 
perse),  each  having  a  characteristic  frequency.  This  effect  is 

1  usually  not  obvious  over  short  distances,  since  the  process 
raiely  has  time  to  produce  complete  frequency  separation. 

b.  Some  frequencies  are  propagated  efficiently  over  rather  long 
distances,  and  others  die  out  quickly.  The  effect  of  this 
phenomenon  is  to  filter  out  some  of  the  original  frequencies, 
leaving  a  signal  characterized  only  by  those  frequencies  that 
are  efficiently  propagated.  In  practical  terms,  the  implica¬ 
tion  is  that  in  some  terrains  and  at  long  detection  distances, 
ail  targets  tend  to  be  characterized  by  wave  trains  exhibiting 
the  same  frequencies.  It  is  obvious  that,  in  such  situations, 
virtually  all  targets  will  look  the  same  insofar  as  the  fre¬ 
quency  composition  of  their  signatures  ie  concerned. 

61.  Variations  in  substrate  characteristics  are  many,  but  the  most  com¬ 
mon  is  stratification.  Nearly  all  areas  are  stratified  to  some  degree,  and 
many  are  divided  very  sharply  into  distinct  layers.  Each  laytjr  tends  to  be 
characterized  by  somewhat  different  propagation  of  seismic  waves,  with  the 
result  that  they  try  to  propagate  independently  through  each.  The  interfer¬ 
ence  caused  by  the  interface  absorbs  energy  and  attenuates  the  wave  rapidly. 

It  should  be  noted  that  stratification  Is  not  a  constant  with  respect  to  time. 
The  most  common  cause  of  a  change  in  stratification  is  a  change  in  moisture 


content  that  can  affect  the  wave  propagation  characteristics  of  the  soil.  In 
addition,  stratification  may  also  change  with  respect  to  small  changes  in 
location,  affecting  wave  propagation  characteristics,  even  drastically. 

Another  cause  of  stratification  changes  is  nonhomogeneous  areas.  These  dis¬ 
parities  occur  quite  often  in  alluvial  deposits  where  a  small  section  of  one 
material  (i.e.,  caliche)  is  surrounded  by  a  different  matrix  material,  pther 
variation  comes  from  faulting  (or  vertical  cracks)  and  changes  in  surface 
geometry  which  provide  more  interfaces  for  complicating  Rayleigh  wave 
propagation. 

.62.  As  the  propagating  wave  reaches  the  geophone,  it  is  carried  on  by 
the  motion  of  the  soil  particles  if  the  geophone  is  properly  placed.  Three 
conditions  of  emplacement  affect  geophone  performance:  (a)  depth  of  emplace¬ 
ment,  (b)  geophone  attitude,  and  (c)  coupling  of  the  geophones  with  the  sur- 
rounding  medium.  Depth  of  emplacement  is  critical  in  measuring  Rayleigh 
waves,  which  usually  have  maximum  signal  close  to  the  surface.  Geophones 
operate  most  efficiently  with  the  axis  parallel  to  the  pull  of  gravity  and 
should  thus  be  leveled  for  use.  Geophones  must  move  with,  the  surrounding 
particles  to  measure  the  motion  and  thus  must  bt  in  solid  contact  with  the 
substrate.  If  not,  the  geophone  will  move  less  than  the  surrounding  parti¬ 
cles,  resulting  in  an  erroneous  measured  signal.  The  preceding  discussion  of 
seismic  wave  propagation  phenomena  is  based  on  an  analysis  by  Benn  and  Smith 
( J  9 75 )  and  more  details  are  available  therein. 

Seismic  Refraction  Data  Correlation 

63.  The'  seismic  refraction  tests  were  conducted  to  define  the  substrata 
of  the  test  site  so  problem  areas  similar  to  those  described  above  are  avoided 
or  at  least  understood.  The  test  site,  because  of  its  location  in  api  alluvial 
valley,  was  expected  to  be  stratified  with  similar  properties  in  each  layer 
and  some  nonhomogeneous  areas.  The  area  was  relatively  homogeneous,  with  sim¬ 
ilar  properties  throughout  the  10-km  test  ar^a,  as  can  be  seen  in  the  seismic 
profile  drawing  of  Figure  10. 

Group  velocities  for  seismic  waves 

64.  The  S— ,  P-,  and  R-wave  group  velocities  calculated  from  the  explo¬ 


sive  and  impact  lists  are  presented  in  the  time  versus  distance  plots  of 
Figure  17  for  comparison  with  the  seismic  refraction  data.  It  should  be  noted 
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Figure  17.  Time  versus  distance  plots  for  explosive-source  tests 

that  the  lower  parts  of  the  curves  are  not  refined  by  extensive  seismic 
refraction  cita  between  0  and  2  km,  and  thus  the  lower  parts  are'  a  biased 
average  of  group  velocity.  The  time  of  arrival  of  each  wave  train  was  cal¬ 
culated  from  the  time  domain  data  of  Plates!  85-87  and  94-99  and  plotted 
against  the  distance  between  the  close-in  geophone  and  the  array.  The  P-wave 
curve  is  similar  to  the  seismic  refraction  data  of  Appendix  B,  indicating  that 
the  P-wave  traveled  in  bedrock  at  depths  over  100  m  with  a  velocity  of  about 
3800  m/sec.  The  S-  and  R-wave  curves  are  shown  only  to  document  group  veloc¬ 
ities  and  their  relationship  to  the  P-wave.  Velocities  from  P-,  S-,  and 

R-waves  are  expressed  by  V  ,  V  and  V  ,  respectively.  These  data  show 

P  5  __  R  • 

that  Vg  is  1277  m/sec  and  VR  is  1364  m/sec  which  gives  V^/Vp*  0.33  and 

V  /V  *  0.94  for  this  site.  The  R-wave  has  a  dominant  frequency  of  3  Hz  in 
K/  b 

the  10-km  test,  indicating  a  wavelength  of  about  455  m,  which  means  that  . 
R-wave  would  be  spread  over  the  upper  ±150  ra.  (This  is  estimated  by  the  fol¬ 
lowing:  dR  -  VR/2f  .)  Thus,  the  long-range  attenuation  data  confirm  and 
support  the  seismic  refraction  data. 
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Seisnic  attenuation 

65.  The  explosive-source  test  data  were  also  used  to  generate  the 
attenuation  curve  shown  in  Figure  18.  This  curve  is  a  plot  of  peak  particle 
velocity  amplitude  (R-wave)  versus  range.  Since  it  uses  data  from  the  1-km 
through  10-km  ranges,  it  is  also  called  a  site  attenuation  curve.  The  curve 
shape  is  another  indicator  that  no  major  anomalies  exist  in  the  substrate. 

66.  Selected  discrete  frequency  vibration  data  were  used  to  prepare  the 
attenuation  curve  presented  in  Figure  19.  This  curve  shows  particle  velocity 
amplitude  versus  range  for  several  different  frequencies.  The  variation  of 
frequency  with  range  can  be  seen  in  the  curves  as  they  flatten  out  as  fre¬ 
quency  decreases.  The  higher  frequency  components  attenuate  faster  than  the 
lower  frequency  components,  which  extend  over  long  ranges.  The  data  were 
taken  at  ranges  ranging  from  5  m  to  1  km  and  relate  quite  well,  showing  that 
site  attenuation  is  relatively  consistent.  This  curve  also  shows  that  the 
data  are  relatively  consistent  over  the  range  of  frequencies,  another  example 
of  the  quality  of  this  data  set. 

Signal  Characteristics 

Mechanical-source  tests 

67.  Vibration.  The  vibration  tests  have  signal  amplitudes  that  cor¬ 
relate  well  over  the  100-m,  500-m,  and  1000-m  test  ranges.  The  only  dif¬ 
ferences  occur  at  the  100-m  range  where  the  amplitudes  vary  between  geophones 
in  the  array.  This  is  probably  because  of  the  significant  difference  between 
the  cylindrical  wave  front  and  a  plane  wave  front  at  100  m  or  because,  for 
continuous  wave  signals,  the  P-,  S-,  and  R-waves  will  add  or  subtract.,  causing 
signal  variations  over  short  distances.  Another  difference  occurs  in  the  tone 
burst  data,  caused  by  variations  in  the  vibrator  output  as  frequency  changes. 
The  data  described  below  are  summarized  in  Table  8. 

68.  The  random  noise  tests  have  signal  characteristics  similar  to  the 
discrete  frequency  vibration  tests.  The  vertical  array  data  from  the  100-m. 
test  have  a  center  frequency  of  25  Hz  and  a  bandwidth  of  10  Hz.  Bandwidth  is 
defined  as  the  width  of  the  band  described  by  the  points  6  db  down  from  the 
peak  amplitude  (one-half  the  peak  amplitude) .  The  500-m  random  noise  test  has 
a  center  frequency  of  15  Hz  and  a  bandwidth  of  15  Hz.  The  - 1-km  random  noise 
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Figure  19.  Attenuation  curves  for  discrete  frequency  vibration  tests 

test  has  a  center  frequency  of  4.6  Hz  and  a  bandwidth  of  1.5  Hz.  The  signal 
appears  to  be  dispersed  to  a  lower  frequency  at  the  i-km  range. 

69.  The  vertical  array  signal  on  the  100-m-frequency  sweep  test  has  a 
center  frequency  of  23.2  Hz  and  a  bandwidth  of  7  Hz.  The  signal  on  the 
500-m-frequency  sweep  test  has  a  center  frequency  of  16.7  Hz  and.  a  bandwidth 
of  1.7  Hz.  The  signal  on  the  1 -kin-frequency  sweep  test  has  a  center  frequency 
of  8  Hz  and  a  bandwidth  of  6  Hz.  This  set  of  tests  shows  the  signal  frequency 
decrease  as  the  range  is  increased,  even  though  the  higher  frequency  signals 
are  input  by  the  source. 

7C.  Impulse .  The  vertical  array  signal  for  the  100-m  impulse  test  has 
a  center  frequency  of  24  Hz  and  a  bandwidth  of  12  Hz.  The  500-m  Impulse  test 
has  a  cfenter  frequency  of  13  Hz  and  a  bandwidth  of  15  Hz.  The  1-km  impulse 
test  produced  a  center  frequency  of  8  Hz  and  a  bandwidth  of  4  Hz.  These  data 
are  listed, in  Table  8. 
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71.  Vehicle.  The  vehicle  tests  produced  signals  of  similar  shape  and 
frequency  content  at  both  speeds.  The  amplitudes  varied  somewhat.  Only  one 
set  of  numbers  is  presented  for  this  discussion.  The  vertical  array  signals 
for  the  100-m  vehicle  tests  have  a  center  frequency  of  12  Hz  and  a  bandwidth 
of  7.2  Hz.  The  500-m  vehicle  tests  have  a  center  frequency  of  2.9  Hz  and  a 
bandwidth  of  2.4  Hz.  The  1-km  vehicle  tests  have  a  center  frequency  of  2.9  Hz 
and  a  bandwidth  of  3.4  Hz.  Vehicle  test  data  are  also  listed  in  Table  8. 
Explosive-source  tests 

72.  Beginning  at  the  1-km  location,  the  explosive-source  test  produced 
a  center  frequency  of  7.2  Hz  and  a  bandwidth  of  3.7  Hz.  This  compares  favor¬ 
ably  with  the  impulse  test  at  1  km.  The  2-km  explosive-source  test  produced  a 
center  frequency  of  3.9  Hz  and  a  bandwidth  of  1.5  Hz  (although  it  is  a  rela¬ 
tively  noisy  record).  The  4-km  test  produced  a  center  frequency  of  3.6  Hz  and 
a  bandwidth  of  1.0  Hz.  The  10-km  test  (number  134)  produced  a  center  fre¬ 
quency  of  3.3  Hz  and  a  bandwidth  of  1.0  Hz,  while  Test  136  (also  10  km,  but 
more  explosive)  had  a  center  frequency  of  3.2  Hz  and  a  bandwidth  of  0.6  Hz. 

The  8-km  test  signal  had  much  wind  noise  and  was  not  included.  It  could 
feasibly  be  used  by  employing  some  frequency  search  filtering  techniques,  but 
this  was  not  done.  These  signal  characteristics  are  also  presented  in  Table  8 
for  comparison  with  other  test  data. 

Background  Noise  Effects 

73.  As  noted  earlier,,  background  noise  was  high  on  certain  days  of 
testing.  Overall,  however,  the  background  noise,  both  seismic  and  acoustic, 
was  relatively  low  except  for  windy  days.  Plates  100-102  illustrate  back¬ 
ground  noise  for  the  three  nondestructive  test  locations,  which  encompass 
three  different  test  days.  The  data  in  Plate  100  for  Test  14  represent  a 
moderate  noise  day  with  engines  off.  The  data  in. Plate  101  for  Test  89  rep¬ 
resent  a  moderate  noise  day  with  engines  running.  The  data  in  Plate  101  for 
Test  132  represent  a  high  background  noise  day.  It  is  interesting  to  note  the 
average  wind  speed  for  the  three  tests:  5.6  m/sec,  5.3  m/sec,  and  10.6  m/sec, 
respectively.  It  appears  that  above  a  wind  level  of  about  8  m/sec,  data 
acquisition  beyond  2  km  is  definitely  impaired  because  the  signal-to-noise 
ratio  may  be  as  low  as  0.1.  Even  with  a  spectrum  analyzer^  it  would  be  dif¬ 
ficult  to  pinpoint  the  signal  of  interest  within  an  extremely  noisy  database, 
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since  the  wind  causes  low-frequency  noise  which  looks  like  actual  data.  How¬ 
ever,  this  problem  could  be  helped  by  deeper  geophone  burial  and  averaging  of 
multiple  sensor  array  data. 

74.  Electrical  noise  is  a  problem  in  some  tests.  The  60-Hz  noise 
introduced  into  the  array  system  produces  problems  when  tests  are  run  above 
50  Hz.  The  signals  can  be  and  were  filtered;  however,  the  filter  does  not 
provide  the  capability  to  separate  real  60-Hz  data  and  the  60-Hz  noise.  The 
electrical  power  also  introduces  harmonics  of  the  60-Hz  noise  in  the  signal 
(especially  120  Hz  and  180  Hz). 

75.  The  background  noise  from  the  generator,  vibrator,  and  impactor 

engines  is  relatively  insignificant  seismically  because  of  the  favorable 

signal-to-noise  ratio  condition.  The  vibratory  or  impulse  data  signal  is  high 

in  relation  to  the  background  noise  caused  by  these  gasoline  engines.  For 

example,  the  seismic  background  noise  measured  5  m  from  the  engines  after 

-4 

Test  89  wqs  about  5  x  10  cm/sec,  while  a  typical  vibratory  signal  was  on  the 

_3 ' 

order  of  5  x  10  cm/sec.  The  background  noise  measured  at  the  array  after 
Test' 89  was  about  1.2  x  10  cm/sec,  while  a  typical  vibratory  signal  at  the 

_5 

500-m  range  was  on  the  order  of  5  x  10  cm/sec.  Thus,  the  signal  appears  to 

be  at  least  an  order  of  magnitude  greater  than  the  noise  caused  by  the 

engines.  The  acoustic  signal  measured  on  the  same  test  at  a  range  of  50  m 

2 

from  the  engines  was  about  10  dynes/cm  ,  while  the  array  microphone  at  a  range 

.  _1  2 
of  500  m  measured  about  1.3  x  10  dynes/cm  . 


76.  Model  predictions  were  made  after  preliminary  site  characterization 
was  completed  but  prior  to  the  initiation  of  testing.  The  model  predictions, 
summarized  in  Figures  20-23,,  are  discussed  individually  below.  -  First,  the 
Rayleigh  wave  phase  velocity  is  plotted  versus  frequency  for  the  preliminary 
site  characteristics  (Figure  20).  Actual  group  velocities  were  about 
660  m/sec  in  the  mid-frequency  range,  compared  to  the  400  m/sec  phase  velocity 
prediction.  The  low-frequency  (3-:Hz)  signal  that  arrived  at  the  10-km  range 
traveled  at  a  rate  of  about  1364  m/sec,  compared  to  820  m/sec  predicted  by  the 
model  since  the  deep  high-velocity  layer  was  not  Included  in  the  model  input. 
Most  seismic  input  by  vehicles  or  other  similar  activity  would  not  have  such 
low-frequency  components.  However,  should  frequency  input  below  3  Hz  become 
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Figure  20.  Rayleigh  wave  phase  velocity  calculated  for  SAT  site 

by  WES  seismic  model 
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Figure  21.  Transmission  coefficient  H  versus 
frequency  predicted  by  model 
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Figure  22.  Transfer  function  T  versus  frequency  predicted  by  model 
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Figure  23.  Transfer  function  T  versus  range  predicted  by  model 


of  interest,  more  extensive  site  characteristics  would  have  to  he  obtained  for 
the  model. 

77.  The  transmission  coefficient  or  R-wave  attenuation  was  predicted  by 
the  computer  model  for  the  test  site,  as  seen  in  Figure  21,  as  plots  of  trans¬ 
mission  coefficient  versus  frequency  for  various  ranges.  This  same  transmis¬ 
sion  coefficient  (or  attenuation)  was  calculated  from  signals  of  the  explosive 
tests  and  is  illustrated  in  Figures  24,  25,  and  26  for  P-waves,  S-vaves,  and 
R-waves.  The  R-wave  values  compare  favorably  with  predictions  but  also  pro¬ 
vide  a  refinement  of  the  model  predictions  in  the  high-frequency  range. 

'  78.  The  site  attenuation  ratio  is  the  relationship  of  the  propagating 
wave  amplitude  with  respect  to  frequency  f^  and  spatial  separation  of  two 
points  (Williams  1981),  given  by 


where 

»  amplitude  associated  with  frequency  f^  at  point  Y,  cm/sec 
■  amplitude  associated  with  frequency  f^  at  point  X,  cm/sec 
Attenuation  ratio  M  is  shown  in  Figures  27-30,  as  calculated  from  the 
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Figure  24.  P-vave  transmission  coefficient  H  versus  frequency  calculated 
from  explosive-source  tests.  Areas  above  about  10  Hr  are  noise  and/or 

filter  contaminated. 


figure  25.  S-wave  transmission  coefficient  H  versus  frequency  calculated 
from  explosive-source  test  data.  Areas  above  about  10  Hz  are  noise  and/or 

filter  contaminated. 
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Figure  26.  R-vave  transmission  coefficient  H  versus  frequency  calculated 
from  explosive-source  test  data.  Areas  above  10  Hr  are  noise  and/or  filter 

contaminated. 
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Figure  27,  Attenuation  ratio  M  versus  frequency  for  100-m  discrete 
frequency  vibration  tests.  Numbers  indicate  location  of  compared  data 


Figure  28.  Attenuation  ratio  M  versus  frequency  for  500-ra  discrete 
frequency  vibration  tests.  Numbers  indicate  locations  of  compared  data 
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Figure  30.  Attenuation  M  versus  frequency  for  l-km,  2-km,  and  4-km 
explosive-source  test  data.  Numbers  indicate  locations  of  compared  data. 


vibration  tests  and  the  explosive  tests  using  the  above  equation.  The  atten¬ 
uation  ratio  M  can  be  calculated  by  the  model  for  an  internal  damping  factor 
of  0.03  by  taking  the  ratio  of  transfer  function  T  for  two  ranges  (from 
Figure  22).  The  values  of  M  from  the  field  tests  were  used  to  calculate  new 
internal  damping  factors  in  the  model  by  changing  the  internal  damping  factor 
for  each  frequency  until  the  proper  value  of  M  was  achieved.  The  results 
are  shown  in  Figures  31-34  plotted  against  frequency,  revealing  several  obser¬ 
vations  that' can  be  made  from  the  data.  The  original  estimated  internal 
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Figure  31.  Internal  damping  versus  frequency  calculated  from  attenuation 
ratio  for  IQO-m  discrete  frequency  vibration  data.  Numbers  indicate 

location  of  compared  data. 

damping  factor  of  0.03  compares  favorably  with  an  average  of  0.02  from  the 
field  data  out  to  1.0  km.  It  should  be  noted  that  at  short  ranges,  a  large 
change  in  internal  damping  factor  produces  a  small  change  in  M  ,  which  means 
that  0.03  produces  M  values  not  unlike  those  from  0.02.  Thus,  future  cal¬ 
culations  for  this  type  of  site  should  include  an  internal  damping  factor  of 
0.02  for  ranges  up  to  l  km.  Both  M  and  the  internal  damping  factor  become 
more  scattered  as  the  range  is  decreased;  for  example,  the  100-ra  and  500-m 
data  curves  are  more  scattered  than  the  1-km,  2-km,  and  4-km  data  curves. 
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Figure  32.  Internal  damping  factor  versus  frequency  calculated  from 
attenuation  ratio  for  500-krn  discrete  frequency  vibration  test  data. 
Numbers  indicate  locations  of  compared  data. 


Figure  33.  Internal  damping  factor  versus  frequency  calculated  from 
attenuation  ratio  for  1-km  discrete  frequency  test  data.  Numbers 
indicate  locations  of  compared  data. 
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Figure  34.  Internal  damping  factor  versus  frequency  calculated  from 
attenuation  ratio  for  explosive-source  test  data.  Numbers  indicate 

locations  of  compared  data. 


This  is  because  at  short  ranges,  small  measurement  errors  in  R-wave  amplitudes 
(or  noisy  signals)  cause  large  changes  in  M  and  larger  changes  in  the 
resultant  internal  damping  factor.  The  2-km  and  4-km  data  indicate  that  the 
internal  damping  factor  may  decrease  even  below  0.02. 

79.  To  weigh 'properly  the  effect  of  frequency  and  range  on  the  computed 
values  of  internal  tjampin£  factors,  it  is  appropriate  to  refer  to  the  proce- 

microseismic  wave  propagation  model  to  normalize,  the 


dure  used  in  the  WEJ 
variables  for  the  tiansmission  coefficient  H 


sion  coefficient  for 
internal  damping  fat 


where 


f 

R 


f  requ 
range 


In  the  model,  the  transmis- 


the  surface  wave  varies  as  a  function  of  kR  and  the 
tor  where  kR  is  defined  as  follows: 


kR 


2  ir  f  R 


ency,  Hz 
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V  =  Raleigh  wave  velocity,  m/sec 

l  v  2tt  -1 

k  =  wave  number  =  ,  m 


The  losses  due  to  the  effects  of  internal  damping  are  in  excess  of  geometric 
losses  due  to  wave  spreading  as  the  waves  propagate  away  from  the  source. 

When  computations  are  made  to  convert  field  measurements  to  values  of  the 
internal  damping  factor,  it  is  appropriate  to  limit  these  calculations  to 
those  data  where  background  noise  and  measurement  errors  are  less  troublesome. 
Thus,  two  limits  must  be  defined:  (a)  the  frequencies  must  be  less  than  those 
defined  as  upper  limits  in  the  data  collection  system  (i.e.,  data  must  have 
good  signal-to-noise  ratios)  and  (b)  the  frequencies  must  be  greater  than 
those  defined  by  a  lower  limit  where  the  internal  damping  factor  is  signifi¬ 
cant.  If  a  value  for  kR  is  assumed  as  30  to  define  this  lower  limit,  a 
frequency-range  table  can  be  assembled  as  follows.' 

Range,  m,  that  must  be  exceeded  for  useful 
Frequency,  Hz  internal  damping  factor  calculation  results 

50 

20 

10 

5 

2 


15 

83 

172 

401 

1815 


This  table  was  used  to  define  the  lower  frequency  limits  for  the  results 
presented  in  Figures  31-34. 

80.  The  site  transfer  function  is  calculated  by  relating  measured  peak 
velocity  amplitude  at  a  particular  range  with  the  load  input  by  the  energy 
source  (Williams  1981).  - 


where 

T^  *  transfer  function  for  frequency 

-  complex  amplitude  of  measured  load  for  frequency  f^ 
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Transfer  function  T  predicted  by  the  model  is  shown  in  Figure  22  plotted 
versus  frequency  for  different  ranges  and  in  Figure  23  plotted  versus  range 
for  different  frequencies.  The  plots  of  T  versus  frequency  for  the  impact 
tests  are  shown  in  Figure  35  and  for  discrete  frequency  vibrations  and  explo¬ 
sive  tests  in  Figure  36.  Figure  35  shows  that  the  consistency  for  the  field 
test  data  is  within  about  5  db  for  the  100-m  geophone  data  from  each  impact 
test  station.  The  curves  compare  favorably  with  the  predictions,  although  the 
values  predicted  by  the  model  fall  slightly  below  the  field  data.  It  should 
be  noted  that  the  site  characteristics  input  into  the  model  did  not  account 
for  the  depth  of  bedrock,  thus  affecting  the  relative  values.  The  change  in 
internal  damping  factor  would  also  have  a  significant  effect.  Background  and 
electrical  noise  (60-Hz)  affected  the  explosive  test  data  above  to  10  Hz,  and 
no  calculation  was  made  for  explosive  source  transfer  function  above  10  Hz. 

81.  Overall,  the  model  did  well  in  making  predictions.  The  basic 
assumptions  used  in  the  model  were  proved  out,  and  the  only  area  In  which 
difference  occurred  was  in  the  low-frequency  realm.  The  change  in  site  char¬ 
acteristics  would  adjust  this  area.  The  slight  adjustment  in  the  internal 
damping  factor  would  make  predictions  better  at  most  frequencies  and  ranges. 
Site  properties  are  quite  important  in  making  model  predictions  and  must  be 
accurate  for  good  predictions  (Lundien  and  Nikodem  1973) . 


100 


Figure  36.  Transfer  function  T  versus  frequency  calculated  for 


PART  VI:  CONCLUSIONS  AND  RECOMMENDATIONS 


Conclusions 

82.  The  set  of  data  obtained  in  this  study  is  relatively  complete, 
encompassing  not  only  close-in  measurements  using  vibratory,  impulsive,  and 
vehicular  sources,  but  also  far-field  impulsive-source  measurements  covering 
signal  variation  in  both  frequency  and  range.  Although  this  analysis  is 
limited,  several  conclusions  are  evident  by  observing  the  data  set. 

83.  Signal  amplitude  at  high  frequency  falls  off  rapidly  with  in¬ 
creasing  range  (paragraphs  48,  50,  64,  66-67,  70,  72,  and  76,  and  Figures  19 
and  26).  This  is  also  illustrated  by  the  fact  that  little  or  no  frequency 
above  50  Hz  was  found  beyond  100  m,  and  the  frequency  sweep  tests  produced 
larger  amplitudes  at  smaller  frequencies  (see  Figure  14).  Also,  only  fre¬ 
quencies  of  about  3  Hz  were  received  at  a  range  of  10  km,  and  in  all  the  ex¬ 
plosive  tests,  only  frequencies  below  10  Hz  were  dominant  in  the  measured 
signal  (Plates  94-98) . even  though  the  explosions  produced  multispectral 
loading. 

84.  Measured  signal  frequencies  are  accurate  within  the  built-in  error 
of  the  FFT  (about  0.5  Hz  in  this  case)  in  comparison  to  the  Input  loading  fre¬ 
quencies  (paragraph  48).  The  frequency  resolution  depends  on  the  length  of 
the  record  and  the  stability  of  the  source.  Tone  bursts  produce  signals  sim¬ 
ilar  to  the  discrete  frequency  vibrations  except  for  vibrator-induced  anom¬ 
alies  such  as  "rounding"  of  the  signal,  the  introduction  of  harmonics,  and  the 
introduction  of  a  large  load  at  the  beginning  and  end  of  the  tone  burst.  In 
cases  where  the  vibrator-induced  anomalies  were  insignificant,  frequency  and 
amplitude  correlate  well  with  the  corresponding  discrete  frequency  test 
results  (paragraph  49). 

85.  Wind  noise  can  be  a  problem  in  making  seismic  measurements.  It 
appears  from  these  tests  that  for  ranges  greater  than  2  km,  wind  sneeds  above 
8  m/sec  cause  a  major  problem  in  signal  measurement.  In  fact,  when  dealing 
with  one  sensor,  it  is  difficult  even  with  a  spectrum  analyzer  to  determine 
where  the  signal  is  when  winds  are  high  and  constant  and  no  precise  method  of 
timing  is  furnished  (paragraph  73).  The  use  of  multiple  sensors,  triaxial 
sensors,  and  a  digital  system  can  help  to  provide  signals  which  can  be  fil¬ 
tered,  cross-sensor  correlated,  velocity  filtered,,  or  otherwise  manipulated 


to  better  define  the  source.  Such  methods  arc  described  by  Aki  (1981),  Allen 
(1982),  Mykkeltveit  et  al.  (1983),  and  Ringdal  and  Huseby  (1982).  However, 
this  type  of  system  would  currently  be  bulky  to  handle  and  complicated  to 
operate  (Benn  and  Smith  1975).. 

86.  Electrical  noise  can  also  be  a  problem  if  a  power  source  introduces 
60-Hz  noise  into  the  signal.  The  60-Hz  noise  can  be  filtered  out,  but  the 
filtering  affects  (he  signal  characteristics  if  the  test  is  run  above  the 
60-Hz  level.  Second  and  third  harmonics  of  the  60-Hz  noise  also  affect  the 
signal  (paragraph  74). 

87.  The  vibration  tests  conducted  at  a  10-percent  offset  frequency  to 
the  3-,  15-,  30-,  and  50-Hz  discrete  frequency  tests  produced  signal  results 
similar  to  the  base  frequency  tests.  The  frequency  difference  is  distinguish¬ 
able  to  the  3-  to  3.3-Hz  tests,  and  the  limitation  on  resolution  appears  to  be 
in  the  FFT  calculation,  although  at  some  point  the  two  frequencies  would  be 
difficult  to  distinguish  (paragraph  48). 

88.  The  vehicle  produced  signals  with  a  low  center  frequency,  and  the 
speed  at  which  it  traveled  made  little  difference  in  signal  amplitude.  The 
vehicle  did  not  produce  a  measurable  signal  beyond  1  km  (paragraph  71). 

89.  The  WES  mathematical  computer  model  can  make  relatively  accurate 
predictions.  These  ttests  have  verified  that  the  concepts  upon  which  the  model 
is  based  are  correct  and  that  the  model's  predictions  are  not  only  qualitative 
but  also  quantitative  in  most  instances.  The  field  tests  verified  the  model 
calculations  of  transfer  functions,  transmission  coefficients,  and  phase 
velocities.  In  most  cases,  the  model  predicted  these  values  within  20  percent 
for  frequency  and  range  (paragraphs  76-78).  Since  the  test  site  was  rela¬ 
tively  homogeneous,  the  model  predictions  are  quite  good  and  sufficiently 
accurate  for  prediction  and  analysis  purposes  (paragraphs  76-81).  It  should 
he  noted  that  quality  predictions  come  from  quality  site  characteristics.  The 
prediction  is  only  as  good  as  the  data  on  which  it  is  based.  The  model  in¬ 
ternal  damping  factor  should  be  refined  somewhat  to  0.02  or  below,  and  the 
site  characteristics  have  been  refined  through  the  complete  testing  (para¬ 
graphs  78  and  80). 

90.  The  seismic  waves  (P,  S,  and  R)  measured  in  these  tests  confirm  the 
data  measured  in  the  seismic  refraction  tests,  and  as  the  range  increased,  so 
did  the  depth  of  travel  and  the  velocity.  This  homogeneous,  deep  overburden 


layer  helps  to  disperse  the  R-wave  into  bedrock  to  a  low-frequency  wave  (par¬ 
agraph  64) . 

91.  It  appears  from  the  lateral  geophone  system  data  that  the  high- 
contrast  interface  of  bedrock  was  so  deep  that  no  enhancement  of  signal 
propagation  occurred  (paragraph  56). 

Recommendations 


92.  Tho  following  recommendations  are  made  as  a  result  of  the  conduct 
of  this  study: 

a.  This  data  set  should  be  used  for  further  analysis  of  seismic 
data.  The  analysis  available  to  date  is  by  no  means  complete 
and  should. be  extended  to  include  pointing  angle  consistency 
and  correlation  with  local  geology,  time  of  arrival  for  mode 
analysis,  cross-correlation  and  coherence  between  wave  types, 
and  extrapolation  of  results  to  other  sites  with  known 
properties. 

b.  The  WES  seismic  model  should  be  used  for  predictions  in  any 
seismic  study  of  surface  wave  propagation.  This  study  affirms 
the  basis  of  the  WES  model  and  its  predictions,  but  also  con¬ 
firms  that  empirical  data  are  valuable  in  verifying  the  model. 
The  internal  damping  factor  should  be  refined  to  0.02  or  below 
for  the  White  Sands  Missile  Range  site,  and  the  site  input 
should  include  a  fourth,  layer  to  simulate  high-velocity  bedrock 
layers  where  appropriate.  The  long-range  data  should  be 
further  evaluated  with  model  calculations  to  determine  whether 
the  internal  damping  factor  of  below  0.02  will  refine  pre¬ 
dictions  for  2-  to  10-km  ranges.  Caution  should  be  exercised 
in  assuming  that  a  particular  site  can  be  predicted  without 
preliminary  data. 

c. .  When  further  long-range  testing  is  done,  a  lateral  geophone 

system  should  be  emplaced,  with  complete  instrumentation  to 
accommodate  signal  analysis. 

d.  Future  tests  should  incorporate  variation  in  both  frequency  and 
range  and  should  include  data  for  as  great  a  range  as  ,is  prac¬ 
ticable. 

e.  This  data  set  is  presented  as  a  generic  study  so  that  the 
results  can  be  applied  to  a  number  of  related  studies.  In 
particular,  these  results  can  be  applied  to  site  security 
analyses  and  activity  detection  studies. 
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Table  1 


Summary  of  Site  Characteristics 


Cone 

Penetrometer  Tests  Relative ' Density 


Station 

Depth 

cm 

Cone  Index 

of  Soil  (Sand) 
g/cc 

Array 

7.5 

100 

1.5  -  1.6 

15.0 

125 

1.5  -  1.6 

12.5 

142 

1.5  -  1.6 

30.0 

250 

1.6 

- 

1.7 

Topography 

37.5 

317 

1.6 

- 

1.7 

0.1  km 

7.5 

100 

1.5 

- 

1.6 

Flat  to  rolling 

15.0 

167 

1.5 

- 

1.6 

22.5 

133 

1.5 

— 

1.6 

Vegetation 

Coverage  % 

30.0 

142  • 

1.5 

- 

1.6 

Sparse/desert 

5 

37.5 

233 

1.6 

- 

1.7 

Grass/creosote 

10 

4.0  km 

7.5 

160 

1.5 

- 

1.6 

Yucca/cactus 

2 

15.0 

165 

1.5 

- 

1.6 

22.5 

160 

1.5 

- 

i  .6 

Weather 

Date 

30.0 

345 

1.6 

- 

1.7 

Snow/ calm 

5-8  April 

37.5 

310 

1.6 

- 

1.7 

Dry / calm 

9-10  April 

10. 0  km 

7.5 

220 

1.6 

- 

1.7 

Dry /windy 

11-13  April 

15.0 

330 

1.6 

- 

1.7 

22.5 

600 

1.7 

- 

1.8 

30.0 

750+ 

i .  7 

- 

1.8 

■ 

37.5 

750+ 

1.7 

- 

1.8 

Station 

Moisture  Content 

Depth;  cm  Value,  7, 

Average  Seismic  Properties 
P-wave 

Depth,  m  Velocity,  m/sac 

0. 1  km 

0-5 

5.4 

0-6 

365-853 

15-30 

9.9 

6-42.7 

853-1525 

l.Okra 

0-15 

6.4 

42.7-91.5 

914-2012 

Table  2 


Summary  of  Mechanical-Source  Tests  Conducted 


Test  No. 

Frequency 

Hz 

Type* 

Date 

Time 

Test  No. 

Frequency 

Hz 

Type* 

Date 

Time 

X 

100- 

-r.  Station 

100-ir.  Station 

fCont . ) 

1 

120 

Discrete 

4/8/83 

1752 

37 

2 

TB  8/8 

4/9/83 

1134 

2 

100 

1757 

38 

i 

TB  8/8 

1138 

3 

50 

180  3 

39 

— 

Random  noise 

1154 

u 

30 

1807 

40 

— 

Random  noise 

1157 

X 

5 

20 

1810 

41 

— 

Random  noise 

1159 

6 

.15 

1813 

42 

— 

Impactor 

1400 

7 

10 

1817 

43 

1405 

8 

7 

1820 

44 

—  . 

1410 

9 

6 

1825 

45 

— 

1415 

X 

10 

5 

1827 

46 

■  — 

1420 

1  1 

A 

1831 

47 

— 

Vehicle  slow 

1514 

12 

3 

* 

1840 

'  48 

•— 

Vehicle  slow 

1518 

13 

') 

1851 

49 

-- 

Vehicte  fast 

1522 

14 

1 

1853  ’ 

50 

-- 

Vehicle  fast 

1524 

X 

13 

120-1 

Sweep 

191 1 

It- 

110 

Offset. 

4/9/83 

0958 

500-m  Station 

1 7 

1 00 

Discrete 

1000 

18 

55 

Of  t  se  t 

1003 

52 

50-1 

Sweep 

4/10/1(3 

1126 

19 

50 

Discrete 

1006 

53 

50- i 

Sweep 

1138 

X 

20 

•16.5 

Discrete 

1010 

54 

50 

Discrete 

1150 

21 

1  5 

Discrete 

1012 

55 

50 

t  B  16/128* 

1  152 

22‘ 

7.3 

01  I  set 

1017 

'*b 

55 

Offset 

1155 

23 

,  3.0 

Discrete 

1019 

57 

30 

Discrete 

- 

1200 

120 

TB  32 

/’  1  28 

10  34 

58 

30 

TB  16/128 

1 202 

T  « 

2  3 

100 

TB  32 

/ :  2« 

1034 

59 

20 

Discrete 

1206 

2b 

50 

TB  18/128 

1045 

b0 

20 

TB  16/128 

1208 

27 

JO 

IB  18/128 

1048 

61 

16.5 

Of  t  set 

1213 

• 

28. 

20 

TB  18/128 

1052' 

62 

’  15 

Discrete 

1215 

>u 

1  5 

TB  18/84 

1 

105  7 

63 

15 

■TB  16/64 

1217 

10 

TB  18/32  . 

1104 

64 

10 

Discrete 

1220  , 

X 

•31 

7 

TB  18/32 

1  108 

65 

■  •  to 

TB  16/32 

1222 

b 

TB  18/32  ' 

1111 

66 

7 

Discrete 

1228 

•  .  *• 

)j  . 

b 

TB  8/32 

1114 

67 

7 

TB  16/32 

1729 

,  34- 

5 

TB  8/18 

1119 

68 

b 

Discrete 

1233 

33 

u 

TB  8/18 

1124 

69 

f» 

TB  16/32 

1217 

i 

38 

3 

TB  8/ 

;  b 

s 

> 

1129 

70 

5 

Discrete 

1242 

’ 

(Continued) 

NOTE:  On  all  vibration,  impart,  and  vchic le  ‘  tests  (except  JiM-129),  data  were  recorded  at  both  the 
,  source  and  the  arrav. 

'*  Test  types  are:  Discrete  -  discrete  frequency;  Sweep  -  frequency  sweep;  Offset  -  offset  fre¬ 
quency  (selected  discrete  frequency  ♦  10%);  Tft  -  tone  bursts  (cycles  on/cycles  off). 


Table  2  (Concluded) 


Test  No. 

Frequent; 

Hz 

Type* 

Date 

Time 

Test  So. 

Krequenc 

Hz 

y 

Tvpe*  Date 

Time 

500-m 

Station 

(Cont . ) 

i 

-km  Station  (Cont.) 

71 

5 

TB  8/ 1  ft 

4/10/83 

1243 

100 

10 

TB  16/32  4/11/83 

1028 

72 

4 

Discrete 

1248 

101 

7 

Discrete 

1033 

73 

4 

TB  8/16 

1249 

102 

7 

Discrete 

1035 

'74 

3.3 

Offset 

1252 

103 

7 

TE  16/32 

1036 

75 

3 

Discrete 

1256 

104 

6 

Discrete 

i0  38 

7  ft 

3 

TB  8/8 

1258 

105 

6 

TB  16/32 

1039 

77 

2 

Discrete 

1301 

106 

5 

Discrete 

1043 

78 

2 

TB  8/8 

1303 

107 

5 

TB  8/16 

1046 

t  'i 

— 

Random  Noise 

1358 

108 

4 

Discrete 

1050 

80 

Random  Noise 

1401 

109 

4 

TB  8/1 

b 

1052 

81 

— 

Vehicle  Slow 

1419  ' 

no 

3.3 

Offset 

11.55 

83 

— 

Vehicle  Slow 

1422 

Ill 

3 

Disc  rete 

1 05  7 

83 

— 

Vehicle  F<,st 

1423 

112 

— 

Random 

Noise 

i  103 

rt4 

— 

Vehicle  Fast 

1425 

1 1 J 

— 

Random 

Noise 

1 107 

85 

— 

Impactor 

1501 

114 

— 

Impulse 

M  12 

8ft 

— 

1502 

115 

— 

1  1  34 

87 

— 

1503 

116 

... 

1135 

88 

— 

1504 

117 

-- 

113ft 

89 

1505 

118 

.  -- 

1  13'7 

119 

— 

Vehicle  Fast 

1153 

1-kra  Station 

120 

— 

Vehicle  Fast 

1154 

' 

121 

— 

Vehicle  Slow 

1155 

90 

30-1 

Sweep 

4/11/83 

0940 

1 22 

— 

Vehlc! 

e  Slow 

1156 

91 

30-1 

Swe  e  p 

0950 

92 

30 

Discrete 

1002 

2- 

-km  iStat ion 

93 

30 

TB  16/128 

1004 

;  94 

.  20 

Discrete 

1008 

123 

6 

•  Di screte 

i  354 

95 

20  ■ 

TB  16/128 

1011 

124  . 

6 

TB 

- 

1356 

96 

'16.5 

Offset 

1020 

125 

5 

Discrete 

, 

1358 

97 

15 

Discrete 

102  3 

126 

5 

TB' 

1  359 

98 

15 

TB  16/64 

1024 

127 

20-1 

■  Sweep 

1402 

99 

10  , 

Discrete 

102  7 

128 

— 

Vehicle  Fast 

1418 

129 

— 

Vehicle  Fast  1 

1420 

*  Test  types  are:  Discrete  -  discrete  frequency;  Sweep  -  frequency  sweep;  Offset  -  offset  fre¬ 
quency  (selected  discrete  frequency  +  10Z);  TB  -  tone  bursts  (cycles  on/cycles  off).1 


j. 


Table  3 


Summary  of 

Explosive-Source  Tests 

Conducted 

• 

2-km  Station 

Test  No. 

Explosive 
weight,  kg 

Explosive 

Charge 

Location,  km 

Date 

Time 

51 

■  12.25 

10 

4/9/83 

1730 

130 

12.25 

10 

4/11/83 

1642 

131 

12.25 

4 

4/11/83 

1824 

132 

'  12.25 

2 

4/12/83 

1432 

133 

12.25 

■  1 

4/12/83 

1602 

134 

12.25 

10  , 

4/13/83 

0859 

135 

12.25 

8 

4/13/83 

1122 

136 

63.50 

10 

4/13/83 

1421' 

NOTE:  On  all  test  shots  (except  51),  data  were  recbrded  at  the  source,  arny, 
and  lateral  systems.  y  -  '  ' 
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AiJ  filtering  was  done  with  a  low-pass  filter. 


Peak  Amplitudes  for  Selected  Impulse  Tests 
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All  filtering  was  done  with  a  low-pass  filter. 

Actual  range  is  1900  m  on  the  8-km  shot. 

10-km  shot  using  63.5  kg  of  explosive.  To  scale  amplitudes  to  an  equivalent  12.25-kg  charge,  multiply 
value  by  0.578. 


LOCAL  i  air  temp  rclative  pyranometer  wind  wind  rain 

^ATE  TIME  ADQVC  SOU.  SlTFACE  HUMIDIFY  ABOVE  SOIL  SURF  SPEED  DIRECTION  GAGE 

JULIAN  CALENDAR  hour  GROUND  TEMP  30  CM  l?.0.  CM  PERCENT  £00  CM  M/SEC  DEG  MM 
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LOCAL  A1P  ICMP  RELATIVE  PYRANOMETER  WIND  WIND  RAIN 

date"  time  adovc  soil  surface  HUMIDITY  aoove  soil  surf  SPEED  direction  cage 

JULIAN  CALENDAR  HOUR  GROUND  TEMP  30  CM  120  CM  PinCENT  200  CM  M/SEC  DEC  MM 
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Table  8 


Signal  Characteristics 


Peak  FFT 

Center 

Band- 

* 

Amplitude 

Frequency 

width* 

Type  of  Test 

Range  m 

cm/ sec 

Hz 

Hz 

Random  noise 

100 

2  x  10_4 

25.0 

10.0 

500 

1  x  10"6 

15.0 

15.0 

• 

1000 

6  x  10~6 

4.6 

1.5 

Frequency  sweep 

100 

4  x  10“3 

23.2 

7.0 

500 

2.2  x  10"4 

16.7 

1.7 

1000 

3  x  10-5 

8.0 

6.0 

Vehicle  (16  km/hr) 

100 

5.6  x  10~5 

12.0 

7.2 

500 

1.8  x  IO-5 

2.9 

2.4 

1000 

*— • 

• 

X 

*— 

o 

1 

Ui 

2.9 

3.4 

Impulse 

100 

1.1  x  10"4 

24.0 

12.0 

, 

500  , 

5.1  x  10*6 

13.0 

1.5 

'  1000 

1.7  x  10“6 

8.0 

4.0 

Explosive  (12.25  kg) 

1000 

2.6  x  10*4 

7.2 

3.7 

(12.25  kg) 

20C0 

2  x  !0~5 

3.9 

.  1.5 

(12.25  kg) 

4000 

1.5  x  10-5 

3.6 

1.0 

(12.25  kg) 

loooo 

4.2  x  10'fe 

3.3 

1.0 

(63.5  kg) 

10000 

4.0  x  10"6 

3.2 

0.6 

Background  noise 

0/5C0** 

1.2  x  10"6 

72.0. 

37.0 

(seismic) 

95/100 

.8  x  10-4 

61.0 

1.5 

495/500** 

r 

o 

X 

in 

48.0 

2.5 

900/1000 

1.2  x  io-4 

10.0 

20. ,0 

■'  800/1000 

-4 

1  x  TO 

5.0  ■ 

15.0 

Background  noise 

■0/500** 

1.3  x  10~l t 

8.0 

14.0 

(acoustic) 

50/100 

7  x  10'1 t 

5.0  . 

11.0 

450/500**. 

io 

5.0 

9.0 

0/1000 

3.8  x  103 

4.0 

6.0 

NOTE:  Background  noise 

test  shows  location  measured/ locat ion  of  test. 

•  0  »  measured  at 

array;  others 

measured  5  m  or 

50  m  from  source. 

*  ,  Width  of  band  described  by  points 

6  db  down  from 

the  peak  amplitude 

(at 

one-half  the  peak  amplitude). 

**  Engines  running  (vibrator,  impactor,  and  generator), 
f  Units  are  dyne:s/cm^. 
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APPENDIX  A:  SEISMIC  ATTENUATION  TEST  PROGRAM 


Introduction 

Background 

1.  Seismic  tests  have  been  conducted  at  sites  all  over  the  continental 
United  States  by  the  US  Army  Engineer  Waterways  Experiment  Station  (WES)  to 
characterize  site  conditions.  Nondestructive  seismic  tests  such  as  those 
which  measure  wave  refraction  are  valuable,  easy-to-run  field  experiments  that 
provide  quantitative  engineering  information  on  surface  and  subsurface  soil 
and  rock.  Experience  has  shown  that  seismic  and  acoustic  signals  are  affected 
by  a  number  of  environmental  source  variables  which  can  cause  large  differ¬ 
ences  in  measured  seismic  response.  However,  there  has  been  little  correla¬ 
tion  of  these  differences,  and  the  effects  of  environmental  constraints  on 
seismic  wave  propagation  cannot  be  predicted  with  confidence  unless  these  ' 
source  variables  are  properly  documented. 

Purpose 

2.  The  purpose  of  this  test  plan  is  to  establish  a  series  of  experi¬ 
ments  that  will: 

a.  Determine  the  reiative  efficiencies  that  seismic  signals  from 
impulsive  and  continuous  wave  sources  propagate  in  natural 
terrain, 

b.  Determine  the  resolution  and  fidelity  of  multiple  frequency 
signals  propagating  in  natural  terrain. 

c.  Supplement  an  existing  database  that  in  turn  will  be  used  to 
validate , theoretical  models  of  seismic  wave  propagation  in 
natural  terrain. 

Approach  '  , 

3.  The  project  will  be  conducted  by  WES  personnel.  Based  on  existing 
data,  a  site  with  appropriate  physical,  seismic,  and  environmental  properties 
will  be  selected.  The  experimental  array  will  be  laid  out  at  the  site,  and 
the  site  will  be  characterized  physically  and  seismically  to  confirm  its 
properties.  An  instrumentation  system  including  geophone,  microphone,  and 
meteorological  sensors  will  be  emplaced  and  monitored  while  seismic  sources 
are  moved  from  station  to  station  along  the  experimental  array.  Measured  data 
will  be  verified  on  site  and  then  processed  and  analyzed  at  WES.  Background 
noise  will  be  monitored  during  the  tests  and  characterized.  The  final  data 


A1 


will  be  used  to  verify  theoretical  model  calculations  and  will  ,be  extrapolated 
to  a  range  of  conditions  through  model  analysis. 

Site  Selection 

Requirements 

4.  The  site  will  be  selected  to  meet  the  general  characteristics  listed 
in  Table  A1  below.  Some  tentative  areas  have  been  located  at  the  White  Sands 
Missile  Range  (WSMR) ,  New  Mexico,  shown  on  Tularos'a  Quadrangle  topographic  map 
distributed  by  the  US  Geological  Survey  (USCS) .  These  areas  provide  com¬ 
mercial  power  and  easy  access,  as  well  as  properties  similar  to  those 
required. 

Table  A l 

Site  Characteristics 


Property 

Description 

Topography 

Flat/rolling 

Vegetation 

Sparse/no  trees 

Weather 

Relatively  dry 

Soil,  surface 

Silt  to  sand  -  dry 

Soil,  subsurface 

Silt/sand  -  cemented 

Bedrock 

Sedimentary  -  D  >  45  m 

Water  table 

D  >  45  m 

Seismic -compress  ion 
wave  velocity 

300  m/sec  (surface)  -  5000  m/sec 
(bedrock) 

Location 

Remote  -  no  obvious  seismic  sources 

Site  layout 

5.  The  site  will  be  laid  out  in  the  array  specified  in  paragraph  10. 

The  actual  location  of  the  site  will,  be  established  by  survey  after  an  on-site 
inspection  and  preliminary  characterization.  The  site  will  be  laid  out  to 
facilitate  ease  of  testing  and  movement  of  equipment.  It  will  be  necessary  to 
drive  to  each  source  location  with  the  vibrators  and  instruments.  Power  will 
be  required  for  the  instrumentation  equipment  at  the  sensor  array,  and  cannot 
be  furnished  by  generator  because  of  acoustic  and  seismic  noise. 


Site  Characterization 


Seismic 

6.  Seismic  refraction  and  vibratory  tests  will  be  made  at  the  array  and 
at  most  source  locations  on  the  test  layout  (Figure  Al).  These  tests  will  in¬ 
clude  measurement  of  compression  and  shear  wave  velocities  with  a  geophone 
package  emplaced  over  a  200-m  length.  This  will  provide  a  seismic  profile 
down  to  not  less  than  50  m.  The  seismic  signals  will  be  measured  with  a 
standard  seismic  refraction  geophone  (see  paragraph  16). 

Environmental 

7.  Each  of  the  seismic  source  locations  will  be  characterized  using  ( 
quantitative  environmental  data  describing  the  substrate  (i.e.,  soils,  under¬ 
lying  rock,  etc.),  topography,  vegetation,  and  meteorology.  The  environmental 
data  will  consist  of  those  factors  that  may  have  significant  effects  on  the 
generation  and  propagation  of  microseismic  waves.  Relevant  environmental 
characteristics  are  as  follows: 

a.  Soil . 

(1)  Soil  moisture  content,  percent. 

3 

(2)  Soil  mass  (or  wet)  density,,  g/cm  . 

(3)  Soil  type  (Unified  Soil  Classification  System  (USCS)). 

(4)  Atterberg  limits. 

(5)  Soil  particle  sizes,  percent  finer  by  weight. 

(6)  Cone  penetration  resistance,  kPa. 

(7)  Plate  loading  test. 

b.  Topography.  Ground  profile  along  the  array-source  line  will  be 

surveyed.  (The  general  topography  will  also  be  documented  by 

USGS  topographic  maps.) 

c.  Vegetation. 

(1)  Type. 

(2)  Height,  cm. 

2  ■ 

(3)  Density,  stems/cm  . 

d.  Meteorology,  (history  and  test  conditions  at  the  array  point). 

(1)  Air  temperature,  deg  C. 

(2)  Wind  speed,  m/s,ec. 

(3)  Wind  direction,  deg. 
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Figure  Al.  Layout  of  sources  and  array 


(4)  Precipitation  (rainfall,  snowfall),  cm. 

(5)  Relative  humidity. 

8.  Soil  samples  for  describing  the  substrates  will  be  obtained  by 
igging  a  pit  (when  possible)  to  a  depth  of  1.0  m  and  hand-augering  a  hole 
pproximately  10  cm  in  diameter  from  the  bottom  of  the  pit  to  a  total  depth  of 

m  (from  the  original  surface).  Soil  samples  from  which  moisture  content  and 
et  density  can  be  determined  will  be  obtained  at  the  surface  and  at  50-cm 
ntervals  to  a  depth  of  2  m.  Bulk  samples  will  be  taken  at  various  depths  (as 
pecified  by  the  field  engineer)  to  obtain  USCS  classification  for  each  iden- 
ifiable  soil  layer.  In. addition  to  the  soil  samples,  cone  penetrometer  read- 
ngs  will  be  obtained  at  the  array  and  at  each  source  location  along  the 
eference  line  at  the  surface  ana  at  7. 5-cm  increments  to  a  depth  of  45  cm. 

9.  Topography  (ground  surface  configuration)  will  be  obtained  by  sur- 
eying  a  detailed  profile  along  the  source-array  line.  The  vegetation  sur- 
ounding  the  site  will  be  described  by  determining  the  plan;  type,  height,  and 
ensity  of  stems  per  unit  area. '  In  addition,  the  dominant  and  codominant 
lants  will  be  determined.  Weather  conditions  will  be  measured  on-site,  and 
ecords  from  the  nearest  official  weather  station  (e.g.,  Holloman  Air  Force 
ase)  will  be  provided. 


Seismic  Data  Collection 


10.  The  tests  will  be  conducted  in  the  array  shown  in  Figure  A1 .  The 
rray  of  calibrated  geophones  capable  of  measuring  1-200  Hz  will  be  emplaced 
t  the  beginning  point  of  the  test  line.  The  array  will  include  five  vertical 
.eophones,  onte  three-dimensional  geophone,  and  one  outdoor  microphone  unit. 

'his  will  be  called  the  array  geophone  system.  Four  calibrated  vertical  geo- 
ihones  will  be  emplaced  near  the  source  to  measure  input  and  close-in  attenua- 
ion.  Also,  a  load  cell  will  measure  input  load  from  the  vibrator,  and  a 
licrophone  will' be  used  to  measure  noise.  These  will  be  called  the  source 
•eophone  system.  In  addition,  a  portable  triaxial  geophone  package  will  be 
unplaced  at  other  locations  on  bedrock  or  on  a  rigid  substratum  to  give  an 
Indication  of  signal  characteristics  at  a  high-contrast  interface.  This  will 
>e  called  the  lateral  geophone  system.  (The  geophones,  microphone,  and 


instrumentation  packages  are  described  further  in  paragraphs  16-18.  Tests 
using  several  seismic  sources  will  be  conducted  at  each  of  eight  ranges  (at 
the  array,  100  m,  500  m,  1  km,  2  km,  4  km,  8  km,  and  10  km).  However,  when  a 
particular  test  no  longer  produces  signals  at  the  array,  that  test  will  be 
omitted  from  the  remaining  ranges.  (The  source  signatures  are  described  in 
paragraphs  13  and  14.)  The  tests  are  listed  and  described  in  Table  A2. 


Table  A2 

Summary  of  Seismic  Test  Sources 


No . 

Description 

Length 

1 

Vibrate  discrete  frequency  sine  waves  of  1,  2,  3,  4,  5,  6, 
7,  10,  15,  20,  30,  50,  100,  120  Hz  at  maximum  force  output 

~3C  sec  each 

test 

2 

Frequency  sweep,  vibrator  at  maximum  force  output, 

1-200  Hz  ' 

3 

Vibrate  discrete  frequency  sine  waves  of  2,  15,  50,  and 

200  Hz  and  vary  frequency  by  1  to  10% 

'  4 

Vibrate  with  tone  bursts  (approximately  8  cycles  of  each 
frequency  in  no.  1) 

-30  sec  w/dead 
time 

5 

Vibrate  with  random  input  (white  .noise)  source;  limited 
bandpass  1-200  Hz 

-30  sec 

6 

Impact  tests:  5  impacts  to  seat  plate  and  calibrate; 
then  dead  time,  followed  by  5  impacts  (explosives  Will 
be  used  as  an  impact  source  if  necessary,  to  obtain  data 
on  the  fax-range  tests) 

-  30  sec  dead 
time  between 
impacts 

7 

M-35  Truck:  normal  -  modified 

-15  sec 

10  mph 

Spectral  array  data 

lli  The  tests  listed  in  Table  A2  will  provide  an  extensive  amount  of 
data  relating  wave  propagation  arid  attenuation  to  input  frequency  and  loading. 
By, varying  the  range  and  the  input  loading,  resulting  data  will  show  the 
effect  of  the  natural  terrain  on  seismic  wave  propagation.  The  vehicle  test 
will  provide  realistic  data  for  comparison.  All  data  will  be  processed  for, 
spatial  information  by  calculating  pointing  angles  for  waves  propagating 
across  the  array.  The  pointing  angle  variations  with  frequency  and  range  will 
be  correlated  with  site  physical  features  (e.g.,  top  graphy,  vertical 


cuts,  subsurface  variations,  etc.),  so  that  similar  results  can  be  extrapo¬ 
lated  to  other  site  conditions.  Note  that  for  spatial  array  analysis,  geo¬ 
phones  from  both  the  fixed  array  and  the  movable  source  will  be  combined  to 
form  aperatures  of  various  geometry  as  range  is  varied. 

Background  noise  data 

12.  The  seismic  motions  resulting  from  both  ambient  and  natural  noises 
and  those  caused  by  the  testing  equipment  will  be  recorded  periodically  during 
and  between  testing.  Specific  tests  to  be  conducted  include  a  24-hr  test, 
recording  data  for  5  min  of  each  hour,  and  a  spatial  array  test  in  which  data 
are  recorded  for  approximately  5  min  at  30  locations  across  the  site.  Thus, 
background  noise  will  be  characterized  both  with  time  and  space  to  allow 
extrapolation  to  other  sites. 

Seismic  Sources 

13.  The  WES  trailer-mounted  shaker  unit  is  a  self-contained  electro- 
hydraulic  vibrator  consisting  of  a  hydraulic  system  and  shaker  assembly  driven 
by  a  gasoline  engine.  The  programmable  function  system  has  a  maximum  force 
output  of  0.91  metric  ton  between  15  and  200  Hz  (output  drops  below  15  and 
above  200  Hz).  The  load  is  distributed  through  a  76.2-cm-diam  base  plate 
which  limits  soil  stress  to  about  20.7  kPa.  Shaker  output  can  be  in  con¬ 
tinuous  sine  waves  at  discrete  frequencies  and  controlled  force  output, 
frequency-sweep  sine  waves  at  controlled  force  output,  pulsed  discrete 
frequencies,  and  random  signals.  The  shaker  was  fabricated  at  WES  using  a 
controller  and  actuator  manufactured  by  Zonic. 

14.  The  impulse  loader  is  a  Dyna  Source  portable  seismic  energy  source 
mide  by  EG&G,  Geometries.  The  system  provides  a  mechanical  energy  source 
powered  by  a  gasoline  engine.  The  system  has  a  piston  weighing  39  kg  that 
drops  1.8  m,  using  a  vacuum  to  prevent  rebound.  The  cycle  time  for  the  system 
is  15-20  sec.  Should  the  impulse  load  be  attenuated  before  reaching  the  array 
on  any  range,  explosives  will  be  substituted  for  the  impulse  load. 

15.  An  M-35  Army  2.5-ton  truck  will  also  be  used  as  a  source  of  seismic 
signals.  It  provides  a  complex  forcing  function  composed  of  both  impulsive 
and  multiple-frequency  components.  Several  runs  will  be  made  at  each  range  on 
a  smooth  road  path  and  a  modified  path  (i.e.,  path  roughness  increased  to  pro¬ 
vide  larger  amplitude  forcing  function). 


Sensors 


16.  The  seismic  refraction  geophones  will  be  Mark  Products  Model  L-10, 
a  standard  vertical  geophone  used  for  seismic  refraction.  The  Model  L-10  has 
a  natural  frequency  of  8.0  Hz. 

17.  The  geophones  to  be  used  in  the  seismic  tests  are  Geospace  HS-10-1 
calibrated  vertical  geophones  with  a  natural  frequency  of  1  Hz,  a  response  of 
1-200  Hz,  damped  at  70  percent  of  critical,  and  a  sensitivity  of  2.95  V/cm/sec 
The  triaxial  (three  mutually  perpendicular  velocities)  geophone  will  be  a  Mark 
Products  Model  L-4C-3D  with  a  natural  frequency  of  1  Hz,  damped  70  percent  of 
critical,  a  response  of  1-200  Hz,  and  a  sensitivity  of  2.35  V/cm/sec.  These 
geophones  will  be  emplaced  according  'to  the  following  specifications: 

a.  Dig  a  hole  approximately  the  same  size  and  shape  as  the  geo¬ 
phone  that  is  to  buried.  The  best  results  will  be  obtained 
with  the  bottom  of  the  hole  flat  and  the  soil  disturbance  kept 
to  a  minimum. 

b.  Place  the  geophone  in  the'  hole  with  the  device  base  in  firm 
contact  with  the  bottom  of  the  hole. 

c.  Take  care  in  backfilling  the  hole.  The  excavated  material 
should  be  replaced  at  as  near  in  situ  conditions  as  possible. 
The  seismic  signal  undergoes  reflection  and  refraction  at  each 
interface.  This  condition  causes  some  of  the  seismic  energy  to 
dissipate.  If  the  soil  is  backfilled  exactly  as  the  soil  in 
situ,  the  number  of  interfaces  between  the  source  and  sensor 
will  be  reduced  by  one.'  Since  this  is  not  altogether  possible, 
the  purpose  is  to  approach  this  homogeneous  state,  which  keeps 
the  properties  on  each  side  of  the  interface  similar  and  the 
energy  decay  to  a,  minimum, 

18.  The  microphone  to  be  used  in  the  array  is  a  Bruel  and  Kjaer  (B&K)' 
Model  4921  outdoor  microphone  unit  with  a  frequency  of  10-20,000  Hz  and  a 
dynamic  range  of  140  db.  This  unit  has  a  built-in  electrostatic  calibrator 
giving  90-db  sound  pressure  level  at  1000  Hz. 


Instrumentation 

Systems 

19.  The  array  geophone  instrumentation  system  will  record  data  at  the 
array.  The  System  will  operate  on  120-V  powet ,  which  cannot  be  provided  by  a 
portable  generator  because  of  the  noise  in  the  array  area.  The  source 
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)hone  system  will  be  portable,  will  always  have  the  same  orientation  to  the 
:ce,  and  will  move  to  each  new  source  location.  This  system  will  be 
:ery  powered  for  convenience.  A  schematic  of  a  basic  instrumentation  sys- 
is  shown  in  Figure  A2.  Both  systems  described  above  are  similar  to  the 
*ing,  except  the  array  system  would  have  a  microphone  added  to  the  input  as 
L  as  a  low-pass  filter  between  the  amplifier  and  the  recorder.  The  lateral 
ihone  system  (paragraph  10)  is  battery  powered  and  has  only  three  channels. 
Ls  a  completely  self-contained  system  with  amplifier  and  analog  tape  re- 
ler.  Its  frequency  response  is  1-200  Hz;  its  natural  frequency  is  1  Hz 
i  a  sensitivity  of  2.35  V/cm/sec. 


Figure  A2.  Block  diagram  of  the  major  components  in  a 
typical  seismic  analog  recording  system 


orders 

20.  Signals  from  the  array  geophone  system  will  be  recorded  on  a 
gamo  Model  3500  14-track  analog  tape  recorder,  recording  at  3-3/4  in. 

5  cm) /sec.  FM-mode,  which  will  give  a  frequency  response  of  DC  to  1250  Hz. 
s  recorder  requires  120  V  AC  power.  The  source  geophone  recorder  will  be  a 
al  Model  STORE  70S,  seven-channel  analog  tape  recorder  operated  in  FM  mode 
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at  3-3/4  in.  (9.5  cm)/sec  (DC  to  1250-Hz  frequency  response).  This  recorder 
will  be  battery  powered,  operating  on  24  V  DC.  . 

Amplifiers 

21.  The  array  system  will  employ  Ithaca  Model  456  amplifiers,  with  a 
bandwidth  of  1-100,000  Hz  and  100-db  gain  range  in  1-db  steps.  Which  operate 
at  120-V  AC.  The  source  system  amplifiers  will  be  WES-made  units,  with  a 
bandwidth  of  DC  to  5000  Hz  and  a  continuous  gain  range  from  0  to  60  db,  which 
operate  at  12-V  DC. 

Filters 

22.  The  array  system  will  use  Krohn-Hita  Model  3323  filters  which  have 
the  option  of  being  low-pass,  high-pass,  or  band-pass.  These  filters  have  a 
frequency  range  of  0.01  to  100, 0Q0  Hz  and  an  attenuation  rate  of  24  db  per 
octave  outside  the  pass-band.  It  is  not  planned  to  filter  the  source  geo¬ 
phones  since  the  signals  at  the  source  will  be  quite  strong  compared  to  back¬ 
ground  noise. 

Calibrators 

23.  Substitution  sine-wave  calibrations  will  be  used  as  an  alternate 
input  to  the  amplifiers  during  system  calibration.  The  sine-wave  calibration 
signal  will  be  monitored  with  a  precision  voltmeter  and  frequency  counter. 

The  sine-wave  calibration  voltage  will  be  compared  to  the  known  geophone 
sensitivities  to. calculate  the  recording  system  sensitivity  in  velocity  units. 
A  time  standard,  IRIG  B,  will  be  recorded  on  both  tape  recorders  for  time  of 
day  and  common  timing  on  both  tape  recorders. 

Spectrum  analyzer 

24.  The  spectrum  analyzer  will  be  used  with  the  array  instrumentation 
system  to  determine  what  and  how  background  noise  data  can  be  filtered  from  ' 
actual  signals.  The  analyzer  to  be  used  is  a  Hewlitt-Packard  3582,  which  hns 
a  frequency  response  of  DC  to  20,000  Hz  and  a  dynamic  range  of  70  ib.  The 
system  also  has  a  plotting  capability  for  permanent  records  made  in  the  field. 

Test  Operations 

Schedule 

25.  A  schedule  of  activities  for  the  period  prior  to  the  field  tests, 
the  actual  field  operations,  and  the  data  analysis  period  is  given  in 
Table  A3. 


Table  A3 
Schedule 


'  Description 

Time  Frame 

Personnel* 

Remarks 

Complete  preparation  of  plan 

21  Jan  83 

EL/GL/ISD 

Approval  of  plan 

Feb  83 

US  Air 

Force 

Approx. 

30  days 

Preliminary  site  survey 

28  Feb-4  Mar  83 

EL/GL 

WSMR 

Site  layout 

14-19  Mar  83 

EL 

Site  characterization  (physical) 

14-19  Mar  83 

EL 

Site  characterization  (seismic) 

28  Mar-2  Apr  83 

GL/EL/ISD 

Testing 

28  Mar- 16  Apr  83 

GL/EL/ISD 

Signal  processing 

18  Apr-30  Apr  83 

ISD 

Model  calculations 

1  Feb-14  May  83 

EL 

Analysis 

18  Apr- 14  May  83 

EL/GL 

Report  draft 

21  May  83 

EL 

*  WES  organizations: 

EL  -  Environmental  Laboratory. 

GL  -  Geotechnical  Laboratory. 

ISD  -  Instrumentation  Services  Division. 

Site  operations 

26.  All  operations  on  site  will  be  done  in  accordance  frith  WES  and  WSMR 
safety  standards.  Contact  will  be  kept  with  Range  Control  at  all  times.  It 
Is  not  anticipated  that  explosives  will  be  used  in  all  tests ,  but  it  may  be 
necessary  to  generate  measurable  signals  in  the  2-  to  10-km  ranges.  If  explo¬ 
sives  are  used,  all  operations  will  be  performed  according  to  an  Explosives 
Safety  Plan  submitted  tp  WSMR,  and  timing  and  firing  will  be  coordinated  with 
Range  Control.  No  damage  is  anticipated  to  the  site  other  than  a  few  small 
holes,  and  these  will  be  filled  upon  completion  of  the  tests. 

Signal  Processing 


Real  time 

27.  One  method  available  to  examine  data  as  they  are  measured  is  an 
oscillograph  record  to  reproduce  time  histories  of  the  recorded  data.  In 
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addition,  the  spectrum  analyzer  provides  a  "quick  look"  processing  capability, 
both  of  real  time  and  recorded  data,  by  displaying  or  a  CRT  screen  amplitude 
siectrum  data  and  phase  spectrum  data.  It  will  also  plot  displayed  data  on  an 
x-y  plotter.  Once  the  data  are  examined,  changes  can  be  made  to  compensate 
for  deficiencies  in  recording  the  source  input. 

Post-test  analysis 

28.  WES/ISD  currently  have  the  capability  of  performing  time  and  fre¬ 
quency  domain  analysis  over  the  entire  length  of  a  test.  Time  history,  auto 
correlation,  cross  correlation,  and  impulse  response  will  be  calculated  and 
plotted  in  the  time  domain.  Fourier  transform  coefficients,  power  spectral 
density,  phase,  cross-power  spectra,  transfer  function,  and  coherence  will  be 
calculated  for  selected  records  and  plotted  in  the  frequency  domain.  These 
frequency  domain  signals  are  acquired  by  time  averaging  Fast  Fourier  Transform 
ensembles  of  512  points  over  any  desired  length  of  a  test.  A  Hanning  window 
will  be  applied  to  frequency  domain  data.  Either  axis  can  be  plotted  linear- 
ally  or  on  1-5  iog  cycles.  The  test  number;  tape  number;  channel  number;  and 
mean,  root  mean  square,  maximum,  and  minimum  values  are  printed  above  each 
plot. 

29.  Utility  software  to  scale  and  reduce  any  test  to  a  new  file 
starting  at  any  selected  sample  is  available.  These  data  will  also  be  plotted 
or  listed  for  examination  prior  to  the  signal  processing  phase.  Utility 
software  will  also  be  used  to  synchronize  data  recorded  on  different  tapes 
that  have  a  common  event  on  each  tape.  The  data  will  be  digitized  at  a  max¬ 
imum  through-put  of  50,000  samples/second,  using  a  12-bit,  ±5-V  ADC. 

Modeling  Analysis 


Description 

30.  The  WES  seismic  model  simulates  the  generation,  coupling,  propaga¬ 
tion,  and  transfer  of  the  seismic  signal  from  the  target  to  the  sensor.  Fig¬ 
ure  A3  illustrates  these  phenomena.  Stress  from  a  target  (vibrator,  moving 
vehicle,  etc.)  applied  over  a  finite  time  causes  the  surface  of  the  ground  to 
deflect.  This  "forcing  function"  coupled  to  the  ground  initiates  particle 
motion  which,  in  turn,  generates  a  seismic  wave.  Seismic  waves  can  propagate, 
in  a  number  of  modes  and  show  a  decrease  in  amplitude  due  to  geometric  attenu¬ 
ation  and  viscous  damping.  The  seismic  wave  amplitude  may  be  further  reduced 
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Figure  A3.  Seismic  wave  propagation 

if  the  wave  travels  near  the  surface  of  a  macrogeometry  feature.  For  example, 
macrogeometry  features  such  as  trenches,  wedges,  or  steps  in  elevation  can 
cause  part  of  the  wave  traveling  along  the  surface  to  be  reflected  in  a  hew 
direction  of  travel  or  converted  to  body  waves,  which  may  or  may  not  travel  in 
the  same  direction  as  the  original  wave.  The  waves  finally  arrive  at  a  mea¬ 
surement  point  (a  geophone)  in  a  pattern  entirely  different  from  their  initial 
pattern.  This  change  is  a  function  of  site  properties  and  range  and  path  of 
propagation. 

31.  To  develop  the  mathematics  to  simulate  the  phenomena  in’  Figure  A3, 
it  is .  convenient  to  consider  the  earth  as  a  homogeneous  medium,  but  such  a 
generalization  does  not  allow  a  prediction  of  seismic  waves  comparable  to 
those  observed  on  seismic  records  measured  in  the  field.  A  more  realistic  and 
useful  assumption  is  one  in  which  the  earth  is  considered  to  be  a  layered 
medium  with  each  layer  considered  homogeneous.  The  approach  taken  in  the  WES 
seismic  model  is  one  that,  attempts  to  develop  the  mathematics  within  the 
framework  of  theoretical  mechanics  to  simulate  the  phenomena  described  above. 
Pretest  predictions 

32.  The  WES  seismic  model  will  be  used  to  predict  the  signals  to  be  . 
measured  by  the  array  geophones.  The  model  will  use  site  characteristics  to 
be  determined  in  the  preliminary  site  characterlzatloa  tests,  with  the  source 
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loadings  described  in  Table  A2  to  predict  the  signals.  These  predictions  will 
be  used  to  set  gains  on  amplifiers  and  other  pretest  calibrations  which  must 
be  done  to  receive  signals  above  system  noise  levels  and  within  the  dynamic 
range  of  the  instruments. 

Post^test  analysis 

33.  The  WES  seismic  model  will  be  used  to  relate  input  loading  to  mea¬ 
sured  output.  This  will  be  accomplished  by  the  following: 

a.  Establish  the  forcing  functions  for  the  input  loads.  For  the 
vibrator  and  Impulse  sources,  the  forcing  function  will  be 
measured  directly  from  the  load  cells.  For  a  complex  vehicle 
source,  the  forcing  function  will  be  computed  using  the  AMC 
mobility  model  vehicle  dynamics  module  (see  Figure  A4  for 
simplified  diagram  of  this  model  for  an  M-35  truck). 

b.  After  establishing  site  characteristics,  the  model  can  predict 
signals  in  both  time  and  frequency  domains  for  comparison  to 
measured  signals  similar  to  those  shown  in  Figures  A5  and  A6. 

c.  One  purpose  of  the  tests  and  model  calculations  is  to  calculate 
a  transmission  coefficient  which,  when  multiplied  by  the  source 
Coupling  coefficient,  will  produce  a  site  transfer  function. 
This  process  can  be  seen  in  Figures  A7-A10.  Once  a  site  trans¬ 
fer  function  is  completed,  signals  for  any  source  can  be 
calculated. 

d.  When  the  calculations  are  complete  for  the  specified  site, 
other  site • conditions  can  be  input  to  the  model  for  extrapola¬ 
tions  of  the  basic  phenomena  to  other  areas. 


Figure  A4.  Schematic  used  to  calculate  forcing  function  for  the  M-35  truck 
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Figure  A6.  M-151  (jeep)  signals,  32  ka/hr,  tine  domain 
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Figure  A9.  Transmission  coefficient 
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APPENDIX  B:  SEISMIC  REFRACTION  STUDY 


Background 

1.  This  study  was  conducted  by  the  Earthquake  Engineering  and  Geo¬ 
physics  Division  (EEGD)  to  determine  the  compression  (P-wave)  velocity  versus 
depth  characteristics  of  a  test  site  located  at  the  White  Sands  Missile 


Range  (WSMR),  New  Mexico.  The  work  was  performed  in  the  period  3-15  April 
1983  by  Messrs.  S.  S.  Cooper  and  S.  0.  Sanders.  In  this  study,  it  was  desired 
to  model  a  typical  alluvial  valley  environment  having  the  following  general 
characteristics: 


a.  Seismic  (P-wave),  profile: 
Depth,  ft* 

0-20 

20-140 

140-250 

250-600 


P-Wave  Velocity,  ft/sec 
1500  to  2400 
2700  to  4300 
5300  i  200 
5800  ±100 


b.  Depth  to  water  table  of  about  150  ft. 

c.  Reasonable  uniformity  in  subsurface  conditions  over  a  planned 
test  distance  of  6.21  miles  (10  km). 

2.  Information  developed  from  earlier  site  selection  studies  had  iden¬ 
tified  certain  locales  within  the  WSMR  as  having  the  general  characteristics 
desired.  The  WSMR  lies  within  the  Mexican  Highland  Section  of  the  Basin  and 
Range  Province, , which  is  characterized  by  a  series  of  tilted  f?ult  blocks 
forming  asymmetric  ridges  or  mountains  and  broad  intervening  basins.  The  area 
of  primary  interest  for  this  study  was  the  Jornado  del  Muerto  Basin  near  the 
Trinity  Site,  where  the  test  line  was  located  as  shown  in  Figure  Bl.  In  gen¬ 
eral,  overburden  in  the  alluvial  fans  and  aprons  closer  to  the  Oscura  Range 
can  be  expected  to  contain  sand  and  rock  fragments  that  grade  coarser  with 
depth,  as  well  as  minor  amounts  of  silts  and  clays.  Near-surface  caliche 
deposits  and  conglomerates  are  also  found  in  those  areas.  Toward  the  center 
of  the  basin,  the  overburden  typically  consists  of  thin-bedded  fine  material 
interbedded  with  occasional  layers  of  coarse  sand.  The  above  conditions  are 


*  To  convert  “eet  to  metres,  multiply  by  0.3048. 


considered  to  be  generally  representative  of  the  typical  alluvial  valley 
environment, 

3.  In  planning  for  the  investigation,  it  was  realized  that  time  and 
funding  constraints  would  not  permit  continuous  refraction  seismic  coverage 
over  the  entire  6.21-mile  (10-km)  length  of  ;he  test  line,  so  it  was  decided 
to  conduct  one  or  more  seismic  surveys  at  each  measurement  location  along  the 
test  line.  Also,  it  was  planned  to  use  field  data  acquired  from  burled  sensor 
arrays  used  in  another  phase  of  this  study,  i.e.,  source  and  main  arrays,  to 
derive  a  field  attenuation  plot  of  peak-to-peak  maximum  particle  velocity 
versus  range.  This  field-derived  plot  would  assist  in  predicting  ground 
motion  response  to  explosive  charges  detonated  at  various  ranges  from  the  main 
sensor  array,  and,  would  also  provide  a  rough  index  of  uniformity  of  ground 
motion  response,  hence  site  conditions,  over  the  length  of  the  test  line. 

Scope 

A.  This  appendix  documents  he  refraction  seismic  and  special  explosive 
tests  conducted,  describes  how  the  various  seismic  refraction  calculations 
were  carried  out,  and  presents  results  of  the  seismic  investigation/site 
characterization. 

Refraction  Seismic  Tests 

5.  It  was  recognized  that  refraction  seismic  survey  lines  up  to  1200  ft 
in  length  would  be  required  to  reach  depths  of  investigation  on  the  order  of 
300-A00  ft.  The  EEGD  has  a  24-channel  SIE  seismic  unit,  with  two  12-channel 
spread  cables  having  50-ft  geophone  spacings;  this  equipment  was  used  in  the 
investigation.  A  100-ft-long  spread  cable  with  10-ft  geophone  spacings 
(12  channels)  was  also  used  to  acquire  near-surf ace  P-wave  data.  A  16-lb 
(7.3-kg)  sledgehammer  and  striker  plate  were  used  as  the  energy  source  for  the 
near-surface  (10-ft  spacing)  refraction  surveys,  but  explosives  were  used  as 
the  source  for  longer  geophone  lines  using  50-ft  spacings.  For  reasons  of 
safety  and  minimum  restrictions  on  transport,  a  two-component  mix  of  Klnepak 
explosive  Was  used  together  with  an  explosive  bridge  wire  (EBW)  detonator  and 
compatible  FS-10  firing  system.  The  Klnepak  components  are  not  classed  as 
explosive  until  mixed  on  site,  and  the  EBW  firing  system  is  insensitive  to 
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accidental  detonation  from  external  influences  such  as  static,  poverline 
radiation,  radios,  etc.  It  had  been  decided  that  relatively  long  line  refrac¬ 
tion  seismic  surveys  vould  be  run  at  the  main  receiver  array  location  and  at 
ranges  of  2.'48  miles  (4  km)  and  6.21  miles  (10  ka).  Near-surface  (10-ft  geo¬ 
phone  spacing)  refraction  surveys  vould  be  used  at  intermediate  measurement 
locations  along  the  test  line.  The  layout  of  the  refraction  survey  lines  at 
the  main  sensor  array,  and  at  the  2. 48-mile  and  6. 21-mile  locations,  is  shown 
in  Figure  B2.  Time  versus  distance  plots  of  refraction  seismic  results  from 
the  various  surveys  conducted  are  vhovn  in  Figures  B3-B10.  Shown  in  these 
plots  are  the  calculated  depths  to  interfaces  between  layers,  the  apparent 
P-wave  velocity  for  each  layer  encountered,  and  the  computed  true  velocity  of 
each  layer.  Also  indicated  on  each  plot  is  the  direction  in  which  the  survey 
was  run,  i.e.,  N-S  indicates  that  the  shot  point  was  at  the  north  end  of  the 
line  and  the  geophone  spread  was  toward  the  south.  Unless  otherwise  noted, 
the  center  of  the  seismic  survey  was  positioned  over  the  survey  stake 
identifying  the  test  line  measurement  location. 

6.  A  refraction  seismic  computer  program,  developed  by  Dr.  H.  N. 

Mooney,  University  of  Missouri,  was  used  to  calculate  the  depth  to  interfaces 
and  the  apparent  and  true  velocities  of  each  refraction  (layer).  Hand  cal¬ 
culations  were  also  performed  to  Validate  the  computer  results,  using  both 
time/distance  and  time  intercept  formulas.  Examples  of  the  quick-check 
hand-calculation  formulas  are  as  follows: 


First  Layer,  D 


Second  Layer,  D^ 

Third  Layer,  D2  •  |  D 

where 

Dg  »  depth  to  layer 

XR  “  distance  along  x  axla  of  plot  to  the  point  where  the  velocity  (slope) 
changes  due  to  presence  of  a  new  refractor 

V  *  P-wave  velocity  of  the  refractors  above  and  below  an  interface 

R 
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Figure  B3.  Main  array,  overburden  teat 
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Figure  B4.  Main  array  (from  previous  geophysical  survey) 
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Figure  BS.  Main  array 
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Figure  B9.  10-ka  position,  H  -  S 
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Figure  BIO.  10>ka  position,  E  •  V 


7.  The  P-wave  velocities  so  calculated  represent  apparent' rather  than 
true  velocities.  To  determine  true  velocity  when  dipping  strata  are  encoun¬ 
tered,  it  is  necessary  to  shoot  forward  and  reverse  seismic  profiles  so  that 
the  dip  and  true  velocity  of  the  layers  can  be  calculated.  The  true  layer 
velocity  is  determined  from  the  expression 

2vyvd 

v  -  - 

T  V  +  V. 
u  d 

where 

-  true  layer  velocity 
Vy  ■  apparent  P-wave  velocity  up-dip 
Vy  ■  apparent  P-wave  velocity  down-dip 

8.  The  above  calculations  are  of  a  simplified  nature;  however,  the  com¬ 
puter  program  used  is  much  more  sophisticated  since  it  considers  all  the  vari¬ 
ables  postulated  by  Snell's  law  of  refraction,  i.e.,  critical  angle  of 
incidence,  velocity  contrast,  layer  dip,  etc.  A  more  detailed  description  of 
these  matters  is  available  from  the  literature  (Herland  1963;  Headquarters, 
Department  of  the  Army  1979) ;  further  discussion  would  be , inappropriate  for 
this  study. 

9.  Also  shown,  in  Figures  B11-B16,  are  P-wave  velocity  profiles  devel¬ 
oped  from  the  rt fraction  seismic  results  presented.  These  plots  depict  the 
layer  dips  calculated  by  the  computer  program  as  well  as  the  computed  true 
velocity  of  each  layer. 
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Figure  Bll.  Main  array  (H  -  S) 
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Figure  B12.  Main  array  (E  -  W) 
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Figure  B16,  10-ka  position  (E  -  W) 
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10.  In  one  phase  of  this  study,  it  was  necessary  to  detonate  explosives 
at  ranges  of  0.62  mile  (1  km),  1.24  miles  (2  km),  2.48  miles  (4  km),  and 
6.21  milep  (10  km)  north  of  the  main  sensor  array.  Ground  motion  data  were 
obtained  at  two  locations  near  the  explosive  source  (at  distances  of  0.06  mile 
(100  m)  and  0.18  mile  (300  m)  from  the  source),  as  well  as  at  the  main  sensor 
array.  These  data  were  checked  in  the  field  and  were  used  to  develop  the 
maximum  peak-to-peak  particle  velocity  attenuation  versus  range  plot  shown  in 
Figure  B17.  From  this  plot  it  is  seen  that  the  attenuation  versus  range  char¬ 
acteristics  were  reasonably  uniform  in  character  and  that  no  severe  aberra¬ 
tions  in  ground  motion  response  were  observed  over  the  length  of  the  test 
line.  Of  course,  much  more  sophisticated  computer  analyses  of  these  data  are 
now  being  performed. 


Summary  of  the  Refraction  Seismic/Site  Characterization  Results 


11.  As  a  preliminary  assessment  of  desired  versus  actual  modeling  of 
site  seismic  characteristics,  one  can  consider  the  composite  P-wave  velocity 
versus  depth  plot  shown  in  Figure  B18.  This  plot  shows  the  desired  seismic 
profile  and  the  average  depth  and  true  seismic  P-wa’*e  velocities  computed  from 
the  field  data  obtained  in  this  study.  The  match  between  desired  and  actual. 
P-wave  velocity  profiles  is  considered  to  be  reasonably  good  for  purposes  of 
this  investigation. 

12.  Another  summary  of  conditions  along  the  test  line  is  shown  in  Fig¬ 
ure  B19,  where  average  layer  depth  and  computed  true  P-wave  velocity  are  shown 
with. respect  to  position  along  the  test  line.  Considerable  license  was  taken 
in  preparing  this  plot  since  continuous  refraction  seismic  coverage  could  not 
be  accomplished,  for  reasons  discussed  earlier.  Water  tables  were  estimated 
from  available  data  (U.  S.  Department  of  the  Interior  1965) ^  The  plot  shown 
is  considered  to  be  in  good  agreement  with  desired  site  conditions  since  the 
velocity  variations  recorded  to  about  180  ft  in  depth  are  relatively  subtle, 
and  local  variations  in  P-wave  velocity  should  not  have  a  significant  influ¬ 
ence  on  shear  or  Rayleigh  waves  that  typically  propagate  at  much  lower  fre¬ 
quency  and  much  longer  wavelengths  over  distances  greater  than  0.6  mile 

(1  km). 
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MAXIMUM  PEAK-TO-PEAK  PARTICLE  VELOCITY,  IN/8EC 


Figure  B17.  Maximum  peak-to-peak  particle  velocity  attenuation 

versus  range 
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Figure  B19.  Seismic  profile  of  test  site 

>  13.  Finally,  the  fact  that  no  major  variations  were  observed  in  the 
plot  of  maximum  peak-to-peak  particle  velocity  attenuation  versus  range,  shown 
in  Figure  B17,  is  evidence  that  subsurface  conditions  are  reasonably  con¬ 
sistent  over  the  6. 21-mile  (10-km)  test  line.  Rad  this  not  been  the  case, 
these  stations  would  have  anomalous  ground  motion  responses  due  to  major 
subsurface  discontinuities  or  features.  It  is  concluded  that  the  site  chosen 
conforms  reasonably  well  with  the  selection  criteria  and  otherwise  meets  the 
objectives  of  the  study.  It  is  also  concluded  that  the  results  obtained 
should  be  applicable  to  other  alluvial  valley  sites  having  similar  site 
conditions. 
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